"ZEITSCHRIFT CH UND MECHANIK 


 INGENIEURWISSENSCHAFTLICHE 2er LIUNGSSEBENEN 


"Band 28 i Es April 1948 | RR Heft 4 


= 


Be2 2 °° 2, HAUPTAUFSÄTZE 
Re) Zusammenhang zwischen ebenen und rotationssymmetrischen 
Grenzschichten in kompressiblen Flüssigkeiten. 

Von W. Mangler in Göttingen. 
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Verlassen zeigt, daß die Berechnung der Geschwindigkeits-, Temperatur- und Dichteverteilung in der 
laminaren Grenzschicht eines symmetrisch ‚angeströmten Rotationskörpers durch eine einfache Transfor- 
motion auf die der entsprcehenden Größen einer ebenen, laminaren Grenzschicht zurückgeführt werden kann. 


‘ Subject of the treatise is ihe computation of the distribution of velocity, temperature, and density in the 
boundary layer of a rotatory body symmetrically floum round by a laminar flow. By a simple transfor- 
mation, this computation cam be reduced to thai of the distribution of ihe corresponding values. appearing 
in @ two dimenstional, boundary layer in the case of a laminar flow. 


ABTOP HOKAsbIBAeT, YTO PAcyeT PacıperereHug CKOPOCTU, TEMNEPaTypE MH INIOTHOCTH B 
HEBHXPEeBOM HOTPAHHYHOM CIIOe CHMMETPHYHO OÖTEKAEMOTO TEN BPAIIEHHA MOIKeT ÖBITb 
CBeNeH IPH IOMOINH IPOCTOTO IIPeoßpasoBaHua K PACYeTy OTUX BONHYRH- Mm cayuan 
INIOCKOTO -HEBHXPEBOTO HOTPAHHYHOTO CIOA. 


I. Einleitung. 


Im folgenden soll über eine allgemeingültige Transformation berichtet werden, mit deren 


Hilfe man die Berechnung der laminaren Grenzschicht an einem Rotationskörper auf die Be- 
rechnung der laminaren Grenzschicht bei der Umströmung eines zylindrischen Körpers zurück- 
führen kann. Die Kontur dieses Körpers, oder was allein wichtig ist, die Druckverteilung längs 
dieser Kontur kann in einfacher Weise aus der Kontur des Rotationskörpers und dessen Druck- 


verteilung hergeleitet werden. Ein solcher einfacher Zusammenhang besteht auch noch, wenn man. 


die Kompressibilität der Flüssigkeit berücksichtigt. 


Dieser Satz ist dem Verfasser seit längerer Zeit bekannt, aber infolge der Zeitumstände 
nieht veröffentlicht worden. Er ist in der Ende des Jahres 1945 in Göttingen erschienenen 
Festschrift zum 60. Geburtstag von Prof. B. Bet z enthalten, ist so aber nur einem beschränkten 
Kreise zugänglich geworden, so daß eine Veröffentlichung an dieser Stelle wohl als angebracht 
erscheint. 

_ Dieser Zusammenhang zwischen rotationssymmetrischen und ebenen Grenzschichten wird 
an dem Beispiel der Überschallströmung gegen eine Kegelspitze erläutert, die bereits früher 
vonHantzscheundWendt[1]auf andere Art behandelt worden ist. Anschließend werden 
mit Hilfe dieser Transformation die inkompressible Grenzschicht an einem gerade angeblasenen 
Kreiskegel und die von der räumlichen Senkenströmung und von der räumlichen Senkenfaden- 
strömung an einer Wand erzeugte Grenzschicht besprochen, wo sich sogenannte ‚ähnliche‘ 
Geschwindigkeitsprofile ergeben. 


I. Die laminare Grenzschicht bei rötationssymmetrischer und bei ebener Strömung. 
Mit. s sei die Bogenlänge längs eines Meridianschnittes durch den umströmten Rotations- 
körper bezeichnet, mit r, der Abstand eines Punktes seiner Oberfläche von der Symmetrieachse, 
so daß die Körperkontur durch die Gleichung r = r, (s) gegeben ist. Die Vorgänge in der lami- 
- naren Grenzschicht an diesem Körper werden bei symmetrischer Umströmung und bei Berück- 
sichtigung der Kompressibilität der Flüssigkeit durch .die Bewegungsgleichung (oder Kraft- 
gleichung), die Kontinuitätsgleichung und den Energiesatz beschrieben. Dazu kommen noch 
die Zustandsgleichung des Gases sowie ein Gesetz über die Abhängigkeit des Zähigkeitskoefli- 
zienten u von der Temperatur T. Mit n sei die Länge des Lotes von einem Punkt der Flüssigkeit 
auf die Wand, mit s, wie schon erwähnt, die Bogenlänge bis zum Fußpunkt des Lotes bezeichnet 
und mit », und v„ die Komponenten der Geschwindigkeit parallel und senkrecht zur Wand. 
Weiterhin sei = c, T die Enthalpie, E = c,T die innere Energie, p der statische Druck, o die 
Dichte und A die Wärmeleitfähigkeit des Gases. Die spezifische Wärme c, und c, sowie die 
Prandtlsche Zahlo = c,u/A sollen als Materialkonstanten betrachtet werden. Dann lauten 
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die oben genannten Gleichungen mit den in der Grenzschichttheorie üblichen Vernachlässi- ; 


gungen: | IE R 
dv dv dp 9a/[ o\ > En 
lu + ng) +). a 
2, (er) + 5, (@ro%) | EN 1 | ed Age % 
RER | di di dp 19/ % Ov,\2 u ”s 
5 ln ln) +) Es oe 


Dabei ist angenommen, daß r, abgesehen vom Bug des Körpers, also abgesehen von der Um- 
gebung des Staupunktes, nicht zu klein ist, daß also der Körper nicht zu schlank ist. Die Zu- 


standsgleichung sagt aus, daß N 
: 2 2 

| DES ER te RE Er (4) 
konstant ist. Der statische Druck p ist im Rahmen der Grenzschichttheorie als von außen ein- 
geprägt anzusehen, also eine gegebene Funktion von s. Man erhält die Gleichungen (1) bis (3), 
indem man die räumlichen Navier-Stokesschen Gleichungen sowie den Energiesatz \ 
(vgl. z.B. Goldstein[2], Ackeret[3]-auf) die angegebenen Koordinaten transformiert 
und dann die üblichen Vernachlässigungen der Grenzschichttheorie anbringt. 

Die Gleichungen für die ebene Strömung unterscheiden sich, wenn man entsprechende 
Bezeichnungen wie oben verwendet, von den Gleichungen (1) bis (3) nur durch das Fehlen des 
Faktors r, in der Kontinuitätsgleichung (2) (vgl. z.B. Goldstein [2]). Man kann also ver- 
suchen, eine Beziehung zwischen den beiden Gleichungssystemen herzustellen. Dazu ist es 
zweckmäßig, noch die Stromfunktion y (s,n) durch die Beziehungen 


oy 0) x 

IHM z, On = — = u RE 
‚einzuführen, da dann die Kontinuitätsgleichung erfüllt ist und nicht weiter beachtet zu werden 

23 braucht. Im folgenden sollen alle Größen, die sich auf die ebene Strömung beziehen, durch 
Überstreichen gekennzeichnet werden. Die Stromfunktion bei ebener Strömung ist dann durch 


en 3 9 
| u = 2: Es Fe 
definiert. Man setzt nun: 
8 
= r2(s BE 
[N ae, Pp(s)=Pp(s), 
h i(s,n)=ils,n), 
role), ss A ri. © 
L ? e(8,Rn)=o(s,n), 
PH, = zvlon), RB,R)=u(sn) 


wobei L eine feste Bezugslänge ist, die z. B. gleich dem Umfang oder der Länge des gegebenen 
Rotationskörpers gewählt werden kann. Die Transformationsformeln für die Dichte und die 
Zähigkeit folgen daraus, daß u allein von ö und o von i und p abhängt. 


Dann gilt 
0 rn 2 TR d.. 0 ee 9 
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was zu ı beweisen Re ei Ber : 


Im Fall einer Seen mit neh ssebane Dichteanderune fehlea in GL. ( 3) BER a « 


das den Druck p enthaltende Glied, das mit der Kompressionsarbeit zusammenhängt. Da o 


\ x konstant ist, kann man dann aus GI. (1) und (2) bezw..aus (9) und (10) allein die Geschwindig- 5 
 keitsverteilung berechnen und nachher -aus Gl. (3) bzw. aus (11) die Temperaturverteilung. Er: 


"Die Transformation bleibt auch in diesem Falle richtig. 
Damit ist der folgende Sat z.bewiesen: | SB 


,. 


: Man kann den Verlauf einer durch die De Ran Pp( s) und die Kontur. r,.(8): 
. gegebenen räumlichen rotationssymmetrischen Grenzschicht dadurch ermitteln, daß man :& SE 


die ebene laminare Grenzschicht zu der ur, 


x 
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berechnet und dann mit Hilfe der Transformation (7) wieder auf die Fanzälichen Größen“ 
umrechnet. 


Durch die Druckverteilung p(s) am Außenrande der Grenzschicht, ist der Verlauf der 


Geschwindigkeit V (s) der Verlauf der Enthalpie :, (s), der Dichte 9, (s) und der Zähigkeil Ä E% 


4a (s) bestimmt. Denn für reibungsfreie Strömungen gilt: 


a dp di, dp» d B 
— Er ne: le re Eee GR 
Eu ‚gs ds 0. ads ds FEB 2 2 


und Us hängt nur von der Enthalpie i, ab (x = c,/%)- 
Bezeichnet, die Zähigkeit und z, die Schubspannung an der Wand n = v so kann man 
den folgenden Aeslonplosen Ausdruck bilden, der nicht mehr von der Reynoldsschen 


° Zahl abhängt: > 5 


VE ee) dl 
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Ein entsprechender Ausdruck, bei dem alle Größen mit einem Querstrich zu versehen sind, gilt 
für die ebene Strömung. Durch die Transformation (7) wird aber aus (13): 


r:(s)ds 
en jr ER 
20. %0 sr: (8) 


Für die Verdrängungsdicke 6, und Impulsdicke , lautet der entsprechende dimensionslose 
Ausdruck: 
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wenn man den immer wiederkehrenden Pe 
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setzt. 


"kleiner. als dr kritische Winkel ist, eine elige 
nach der wieder Potentialströmung herrscht, 
dureh die Kegelspitze gehenden Geraden ist der Are 


stant. Man hat also in Gleichung (1) bis (3). en 
proportional ist. Man kann demnach nach bis’ (11) a aminare G! 
ebene Platte bei demselben Druck ?= 9, deneiken: Dichte %= = 0a ar 
!%a = Hi, berechnen, was bereits von Hantzsche und Wenc dt [5] « 
die Ergebnisse auf den Kegel übertragen. Der Umrechnungsfaktor für i 
- lautet nach GI. (14) bei demselben ar, XOB.B Und BEE Dose 
A Rn 1a 
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7 
so daß die abet am Kegel Bene der Platte vergrößert ist, 


e- ® 
Grenzschichtdicke am Kegel nach Gl. (16) gegenüber der Platte um den Faktor 1//3 ve ei 
ist, ein Ergebnis, das auf andere Art schon von Hantzsche und Wend +1 1 gefunden. 
wurde. ; 


IV. Dir Kreiskegel hei Kkerenker Strömung. RER, E 
Als weitere Anwendung der Transformation (7) soll noch die lJaminare Grenzschiht n 
einem Kreiskegel in einer inkompressiblen Flüssigkeit bei symmetrischer Anblasung behandelt 
- -  werden.. Die Geschwindigkeit V am Be | E 
ve Außenrande der Grenzschicht ist 
einer Potenz der Bogenlänge s längs 
der Kegelerzeugenden proportional: 


F nd a = 


Der Zusammenhang des Exponenten Be 
m mit dem halben Ditnungseiakel % ER 
des Kegels („= s-sin d%) ist aus a 
Abb. 1 zu ersehen). Die zugehörige < 
ebene Potentialströmung ist wegen (7) 


gegeben. Das ist eine ebene Strömung 
gegen einen Keil vom Öffnungswinkel 5 
Br, wobei 


A re 


ist, wie man sich mit Hilfe der konfor- 
men Abbildung des Keiles auf eine 
3, 0° gerade Wand (8= 1) leicht überlegen 


0° 30° 60° 90° 120° 
— o 


x ‘kann (Abb. 1, vgl. auch [6)]). 


Nun führt aber die Berechnung 
[v) 0.25 ) 17. 
ER % 5 der Grenzschicht an einem Keil auf die 
7 


folgende gewöhnliche Differential- 
Bild 1. Vergleich der räumlichen Strömung gegen einen Kegel gleichung 


mit der ebenen Strömung gegen einen Keil. a +a EL +ß (1 ae, 1a) es) (22) 


Bern ‘) Verf. verdankt dieses Bild Herım R. Leuteri tz, der die symmetrische Potentialströmung gegen eine solche 


Kegelfläche mit Hilfe von Kugelfunktionen berechnet hat. . 


gl. z.B. [6]) | 


"Dabei bedeuten & und N dimensionslose Koordinaten die mit der Str | ER. 
> Or P5 Brak und ’ omfunk Di, 
' und dem Wandabstand in folgender Weise zusammenhängen: Br. 
ö Zn ER TER 2 Sg j 


new... ee 
2 ER er en | und | a a er =) R >eE 
0. ist. Daes auf einen gemeinsamen positiven Faktor von & und ß nicht ankommt und der Rei- 3 £ 
Be: winkel zwischen 0 und liegt, darf hier «= 1 und m = ß/ (2— ß) gewählt werden. Be: 
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Bild 2. Laminare Geschwindigkeitsprofile bei beschleunigter Strömung. 
Die Lösungen der Differentialgleichung (22) mit den Randbedingungen (23) sind von 


‚Hartree [7] tabuliert worden. Aus ihnen erhält man mit Hilfe von Gl. (7) sofort die Ge- 
schwindigkeitsprofile am Kegel. Die Wandschubspannung und die Grenzschichtdicken findet 


man aus: 
# er; z dr. N n d8-N 
E ea Er e He Fed eh BRATEN 4* = ' 1 - EN dn—= —— 
e en u UT ESr [EN an nz ©). 


wobei z;, 6; und ö, durch (13) und (15) definiert sind und der Umrechnungsfaktor nach (17) 
den Wert N = /3 bekommt. 

Insbesondere wird beim gewöhnlichen rotationssymmetrischen Staupunkt d, = = also 
m—1,m= = a=1, ß=0,5. Das vn Homann [8] angegebene rotationssymmetrische 


Staupunktsprofil enthält leider einige numerische Ungenauigkeiten. Man hat in diesem Fall 

ebenso wie bei der ebenen Staupunktsströmung (« = 1, ß = 1) gleichzeitig eine strenge Lösung = 
‘der Navier-Stokesschen Gleichungen. Beide Geschwindigkeitsprofile sind zusammen 7 
mit einigen Zahlenwerten zur Berechnung der Wandschubspannung und Grenzschichtdicke in 
Bild 2 dargestellt. Es sei noch erwähnt, daß die zur räumlichen Staupunktsströmung gehörige 
ebene Ersatzströmung (m Be z an einem Keil mit dem Öffnungswinkel Br = 5 auf- 
tritt. 

Für die Übertragung der Keilströmung auf die Kegelströmung besteht nur die Einschrän- 
kung, daß der Öffnungswinkel 9, des Kegels nicht zu klein sein darf, damit die in der Grenz- 
schichtgleichung (1) enthaltene Voraussetzung erfüllt ist, daß r} - Vu/Lv groß gegen 1 ist. Außer- 
dem ist ebenso wie bei der Keilströmung ein gewisses Gebiet um die Spitze s= 0 herum von 
der Betrachtung auszuschließen. 
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V. Die räumliche Senke und ee a Senkentaden. 
Weiterhin soll die laminare Grenzschicht am He | | 
den, in dessen Mitte sich ein kleines Loch befindet, so daß die Flüssigkeit en 
das Gefäß hinreichend groß und die Öffnung hinreichend klein, so kann man die 
Loches abgesehen von der Umgebung der seitlichen Wände und des Loches durch die 
einer räumlichen Senke ersetzen. Ohne Reibung hat man dann am Gefäßboden eine 
digkeitsverteilung der Form 


BE ner BE en 


Da dsstets in Strömungsrichtung positiv gewählt worden ist, ist #8 —r, in diesem Falle 
negativ. Nach (7) hat die zugehörige ebene Geschwindigkeitsverteilung die Form 5 
2 > rn 

v ” ie nA 2 28 5 

V = TE ’ nn 5 in RE ( ), - 


so daß ß/x = —4 ist und wegen (26) = — 1,8 4 zu wählen ist, Das Geschwindigkeils- 
profil ist zusammen mit einigen Zahlenwerten in Bild 2 zu finden. Die Auftragung über 


nYa=nyı DL ARBE i ar 


v ds » ds 


ist für Grenzschichten mit Druckabfall besonders geeignet, da sich die Profile für a = +1 und 
& = — 1 hier gemeinsam ohne Überschneidungen bequem darstellen lassen, 

Bei der entsprechenden zweidimensionalen Aufgabe, die Grenzschicht in einem Gefäß zu 
berechnen, wo die Flüssigkeit durch einen unendlich langen, schmalen Spalt abläuft (ebene 
Senke), ist die Geschwindigkeit Y zu (—$)"! proportinal, so daß a= 0 und = list, Das 
Geschwindigkeitsprofil (Bild 2) wurde bereits vonK.Pohlhausen [9] im Jahre 1920 an- 
gegeben. ia 

Schließlich soll noch die Grenzschicht an einer Wand betrachtet werden, bei der die Strö- 
mung außerhalb der Reibungsschicht von einem Senkenfaden erzeugt wird, der auf der Wand 
senkrecht steht. Unter einem Senkenfaden soll eine mit Senken konstanter Stärke belegte Ge- 
rade verstanden werden. Dann ist die Geschwindigkeit V zu 9°" = (—s)"! proportional, also 
= —1l,ß= 1 zu setzen. Das Geschwindigkeitsprofil ist ebenfalls in Bild 2 dargestellt®), Aus 
dem Bild ist das verschiedenartige Verhalten der Lösungen der Differentialgleichung (22) für 
großen Wandabstand gut zu erkennen, Wie schon in [6] gezeigt worden ist, gilt die asympto» 
tische Entwicklung (28 +a >): 

EIN) 
= (/»6 i m -} 4) art füra = +1 | 
| Nu nY® are | > 
% Er 


On +4) für a» 


(28). 


In allen diesen Fällen ist die Lösung der Differentialgleichung (22) durch die Randbedin- 
gungen (23) eindeutig bestimmt, wenn man noch verlangt, daß die Grenzschichtdieke endlich 
sein soll (vgl. [6]). Bei der durch Umkehrung der Senkenladenströmung entstehenden Quell- 
fadenströmung ist x = 1, ß = —1 zu setzen. Mills[10] fand in diesem Fall unendlich viele 
Lösungen. Alle diese Geschwindigkeitsprolile haben „Übergeschwindigkeiten“, d.h. die Ge- 
schwindigkeit im Inneren kann größere Werte annehmen als am Außenrande der Grenzschicht, 
Dasselbe gilt für die räumliche Quellströmung (a = 1, = -—A), die Umkehrung der räum- 
lichen Senkenströmung. Erwähnt sei der Vollständigkeit halber noch, daß im Fall der ebenen 
Quellströmung (x = 0, ß = I) überhaupt keine „ähnlichen“ Lösungen der Grenzschicht- 
gleichungen existieren. 

Abschließend sind in Bild 3 einige Kurven dargestellt, die die Abhängigkeit der Wand- 
I a a ee AR 
schubspannung und der Grenzschichtdicken von dem Parameter AA) = = 1, eigen, Der 
„Formparameter“ B/a ist ebenfalls angegeben, Um Unendlichkeitsstellen zu vermeiden, ist eine 
nichtlineare Skala verwendet, so daß die Skala lür linear wird, 
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°) Die Lösung wurde dem Vorf, von Horn K, Wioghardt In dankonsworter Welse zur Vorfügung gostollt, 
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Das Vorzeichen von « und ß ist so zu wählen, daß (2%-—ß)  »/L positiv wird, wobei zu beachten. 
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Bild 3. Abhängigkeit der Wandschubspannung und der Grenzschichtdicke vom 
Druckgradienten bei „bhnlichen‘‘ laminaren Geschwindigkeitsprofilen. , 


ist, daß ds stets in Strömungsriehtung positiv ist. Bei reibenden Strömungen kann man ja 
im Gegensatz zu reibungslosen Strömungen die Richtung nicht umkehren. 


VI, Zusammenfassung. 


Über den Zusammenhang von drehsymmetrisch räumlichen und zweidimensionalen Grenz- 
schichten ist der folgende Satz bewiesen und seine Anwendung an Beispielen erläutert worden: 

Die Berechnung der Geschwindigkeitsverteilung und der Temperatur- und Dichtevertei- 
lung in der laminaren Grenzschicht an einem symmetrisch angeblasenen Rotationskörper kann 
durch eine einfache Transformation (7) auf die Berechnung der entsprechenden Größen einer 
ebenen, laminaren Grenzschicht zurückgeführt werden. Beide Vorgänge lassen sich mathema- 
tisch durch dasselbe Gleichungssystem mit denselben Bandbedingungen beschreiben. 


YIL Sehrifttum., 


[1] W Hantzsche u. H Wendt: Die laminare Grenzschieht bei einem mit Überschallgeschwindigkeit 
angeströmten nichtangestellten Kreiskegel, Jahrbuch 1941 d, Dt. Luftfahrtforschung, 8. I 76. 
8, Goldstein and others», Modern development» in fluid dynamies. Oxford (1938). 
J. Aekeret: Gasdynamik in Handbuch der Physik, Bd.7, Berlin (1927), 
# A, Busemann, Drücke auf kegelförmige Spitzen bei Bewegung mit Überschallgeschwindigkeit. Z. 
angew. Math. Mech, Bd. 9 (1929), B. 496. , 
[5] W.Hantzscheu. H, Wendt: Die laminare Örenzschicht der ebenen Platte mit und ohne Wärme- 
über; unter Berücksichtigung der Kompressibilität. Jahrb, 1942 d. Dt. Luftfahrtforschung, S.I 40. 
[6] W. Mangler: Die "ähnlichen. Lösungen der Prandtlschen Grenzschichtgleichungen. Z. angew. Math. 
Mech. Bd. 23 (1943), 8, 241, 
[7] D,B, Hartree, On an equation oceurring in Falkner and Skan’s approximate treatment of the equa- 
tions of the boundary layer, Proc, Cambridge Philos, Boc, Bd, 33 11 (1937), 8. 223. 
[8] #,. Homann: Der Einfluß großer Zähigkeit bei der Strömung um den Zylinder und um die Kugel. 
Z. angew. Math. Mech. Bd. 16 (1936), 8, 153, 
(9) K.Pohlhaus en: Zur annäherungsweisen Integration der Differentialgleichung der laminaren Grenz- 
schicht, Z, angew. Math Mech, Bd.1 (1921), 8. 252, 
[10] B. H. Mills, A note on some boundary layer veloeity profiles, Journ. Aeron. Bei. Bd. 5 (1938), Nr. 8, 
8, 325, 


+ 


Hingegangen ; 39. Beykemiber 1947. 


A 4 


r . 


Von Alexander Kromm, Berlin. | 

Es wird die Ausbreitung der elastischen Wellen in ebenen Scheiben mit kreisrundem Loch untersucht, 
die Be plötzlich aufgebrachte drehsymmetrische Störung am Lochrand ausgelöst we a 
binierte Anfangs- und Randwertproblem wird mit Hilfe der Laplace- Transformation behandelt. 
der Übersetzung der Lösung im Unterbereich in den Oberbereich ergibt sich eine Integralgleichung erster 
Art vom Volterratyp, die sich leicht numerisch auflösen läßt. h a r 

‚ion of elastic waves in plane discs provided with a circular hole is investigated in 
en a by a sudden Fe mmetric disturbance at the edge of the hole. This 
of both initial and boundary conditions is treated by aid of the Laplace transformation. From t Rast 
formation of the solution of the inferior range into the Ba there follows an integral equation of 
first kind of the Volterra type that can be easily numerically. 

B »To& pa6ore necsenyerca PacıpocTpaHeHHe YIPyTHX BONH B INIACTUHKAX C KPYTOBEIM 
OTBEPCTHEM, KOTOPEI® BEISEBAITCH BHESAIHEIM TIOABJICHHEM BOSMYINCHHÄ, PABHOMEPHO Pac- 
mpeneieHHEIX IIO KPAsM OTBePCTHA. ITA Banaya C HAYAJNBHEIMH H KPaeBEMH 
pemaerca npm nomomu ıpeoßpasopanun Jlarmaca. Ilpu oöpamennu pemenns oToßpameH- 


HOrO NuUhepeHNMaNbHOrO VPABHeHHA B O0NACTL UCKOMOK PYHEIHH TONyYaeTcA HHTEerPaNBHoe 
ypasuenne mep5oro pona runa Boupreppa, Nerko pemaımeeca UMCHeHHEIMH MeTONaMH. 


I. Einleitung. Aufgabenstellung. 

In letzter Zeit haben die Probleme der Wellenausbreitung in steigendem Maße das Inter- 
esse der Elastiker gefunden !). Als ein Beitrag zur weiteren Kenntnis dieser Erscheinungen 
soll hier ein solcher Ausbreitungsvorgang in zwei Dimensionen untersucht werden, und zwar 
die Ausbreitung von Störungen in einer unendlich ausgedehnten ebenen Scheibe mit einem 
kreisförmigen Loch vom Halbmesser r = a, an dessen Rand drehsymmetrische Kräfte plötzlich 
aufgebracht werden. Wir beginnen mit der Aufstellung der Differentialgleichung dieses Vor- 
gangs: 

Die Bewegungsgleichung für eine radiale Verschiebung eines Scheibenelements der Massen- 
dichte o heißt nach Abb. 1 7 

f) u 
a rap dr =, dpAr—rapdre —: a a (1). 
Die darin vorkommenden Spannungen o,, o, lassen sich nach dem Hookeschen Gesetz für 
einen drehsymmetrischen Vorgang durch die radiale Verschiebung % ausdrücken: 
‚[9u U I “) 
au An 
wobei wir die Abkürzung 


HR WEL ee (3) 


eingeführt haben (#E=Elastizitätsmodul, » = Querkontraktionszahl). Setzt man das Elastizi- 


tätsgesetz in die Gleichgewichtsbedingung ein, so ergibt sich eine partielle Differentialglei- 
chung für «: 


92 u 1 ou u 1 u : 
ee ee er BT EEE - 
wobei 


ei VE’/o a Da ee (5) 
die Dimension einer Geschwindigkeit hat. 
Wir suchen Lösungen der Differentialgleichung (4), die den Anfangsbedingungen 


V=l; ef) =Dfür = 0 undalle >, a \ 


genügen. Ferner haben wir am Lochrand r = @ eine Vorschrift für alle Zeiten t > 0 zu machen. 
Dabei wollen wir zwei Fälle unterscheiden, die zur Eingrenzung eines praktischen Problems 
nützlich sind: Entweder es wirke von t= 0 an auf den Lochrand plötzlich ein radial nach außen 
gerichteter Druck: 

ou 27 


Re = 4) für ra gandt 202 Sum (7.1), 


nn S. z.B. W. K ucha L8 ki, Unstetigkeitsstellen in einem bewegten Kontinuum. Z. angew. Math. Mech. 21 (1941), 
S. 152. — W.# lügge ‚ Die Ausbreitung der Biegungswellen in Stäben. Z. angew. Math. Mech. 22 (1942), S. 812. — 
G. Sonntag, Biegestoß gegen eine Kreisplatte. Forschung 14 (1943), S. 137, 
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Zur Ausbreitung von Stoßwellen in Kreislochscheiben. 


TE AR da PT 


INT 


alas Bi 


T Den hi aa ai ae r abe e Ko 7 MT 


Pe a 


z an 


a 2 N, = ve =. für »= a und 20 ET rd & 0 


sen dt 


An Stelle. einer zweiten Randbedingung tritt in beiden Fällen die Forderung, daßufürrr>o 
2 eine reguläre Funktion sein soll. = Bi 
es In der vorliegenden Veröffentlichung werden wir nur den ersten Fall behandeln. Fir den 
Ser zweiten Fall, der sich nach derselben Methode wie der erste lösen läßt, werden wir kurz die Er- 


ne ei En 


gebnisse in der Zusammenfassung erwähnen. 


II. ehe - 


- Wir Ep ein Problem zu lösen, bei dem hinsichtlich t gewisse Anfangsbedingungen (s. 
Gl. (6)), hinsichtlich r gewisse Randbedingungen (s. Gl. (7.1) oder (7.2) und die Regularitäts-- 
. bedingung bei r = oo) vorgeschrieben sind. Eine geeignete Methode zur Behandlung solcher 


‚Probleme beruht an der Anwendung der Laplace- Transformation auf die Zeitvariable t: 


| Fo=jemrga| er 
- oder abkürzend . Fo)-—>f) | 
mit 32 Ir ERROR LEERE SA ter 


"Man Kind 5 (t) die Oberfunktion und F (?) deren Unterfunktion?). Wir bezeichnen die Unter- 
. funktion von 4 (r, £) mit U (r, p) wg 
} Ur) > ul) ne ER TE (10); 
Von großer Bedeutung bei der Anwendung der Laplace-Transformation ist der Ernten, 
daß der Differentiation im Oberbereich eine Mulitplikation mit der Variablen p im Unterbereich 


| £ : 2 .; a 
entspricht und daß in der Laplace-Transformierten einer Ableitung De die Anfangswerte der 


Funktion %(r,t) und ihrer niederen Ableitungen auftreten ®): 


vu °u (r, 0) „lu (r,0) 

ah v-1 Be a hr ee 
TTZ +:-P U(r, p) — pr!u(r, 0) —p FR ...— p°. apa (Eh) 
Damit wird die partielle Differentialgleichung (4) in unserem Falle in eine gewöhnliche über- 
gehen und außerdem werden die Anfangswerte in die transformierte Differentialgleichung selbst 
eingehen, so daß sie dort automatisch berücksichtigt werden. 

Für die Transformierte der zweiten Ableitung nach der Zeit, die in der Differentialglei- 
chung (4) vorkommt, erhalten wir bei unserer Aufgabe wegen der Anfangsbedingungen (6) 
einfach - . - 


ÜBEN INT ER NEN ee er} 

und durch die Laplace- Transformation der Differentialgleichung (4) erhält man für U(r, p) 
tatsächlich die gewöhnliche Differentialgleichung | 
@U 10U iz ak 

er Tr dr \a'n 


Sie unterscheidet sich von der Besselschen Differentialgleichung nur in dem Koeffizienten 


\0=02 2 REN (13). 


von U. Ihre Lösung könnte man in Besselschen Funktionen vom Argument ier auf- 


schreiben ; einfacher ist es jedoch, die Lösung in reeller Form unter Benutzung der Watson- 
schen Funktionen *) K,und I, darzustellen, die im folgenden Zusammenhang mit den Bessel- 
schen und Hankelschen Funktionen stehen: 


1 ; —lysi 
Kl) Zmie”” Hin), Peer 


2) In der mathematischen Literatur wird gewöhnlich die Oberfunktion mit großen, die Unterfunktion mit kleinen Buch- 


staben bezeichnet, vgl. G. Doetsch, Theorie und Anwendung der Laplace-Transformation, Berlin 1937 (im folgen- 
den kurz zitiert als Doetsch mit Seitenzahl). Wir müssen hier von dieser Schreibweise abweichen, um die für die Ober- 
funktionen gebräuchlichen Bezeichnungen der Mechanik beibehalten zu können. 

®) vgl. Doetsch, 8.153. 

48.G.N. Watson, A treatise on the theory of Bessel functions, Cambridge 1922, S. 77 und 78 (im folgenden kurz 
zitiert als Watson mit Seitenzahl). 


ode es wird den Runlpunkten von = = 0 3 eine radiale Bewegung mit konstanter Geschwin- 
BT Bet 0: Bun % 


ITORR | | ee 2% 
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Mit K, und /, erhält man die Lösung von (13) in der Form: 


wobei die Integrationskonstanten noch von p abhängen können. Sie werden aus der Rand- 


bedingung (7.1) oder (7.2) (die in den Unterbereich übersetzt werden muß) und aus der 


Regularitätsbedingung für —oo bestimmt. Wir gehen in diesem Abschnitt nur auf letztere 


ein, da sie den beiden zu untersuchenden Fällen gemeinsam ist. Zu diesem Zweck schreiben wir 


die asymptotischen Formeln für K, und J, für r> @ auf): 


il 2: 
x, (?) = RGn c 


C 2pr 
pr pr , 3ni 
PEN LE RE N. 
A ) 2. a 45 


Da nach (9) Rep>0 sein soll, so geht für r>oo auch I, —, und aus der Regularitäts- 


bedingung für U folgt, daß im Lösungsansatz (14) B=0 sein muß. Die Lösung von U 
enthält somit nur X, (pr/e): 


U(r,p)=A(p)K, 1 ee (15). 


Es bleibt noch A(p) aus der Randbedingung bei r = a zu bestimmen und die gefundene Lösung 
für U in den Oberbereich zu übersetzen. 


III. Lösung der Aufgabe, wenn die Spannung am Lochrand gegeben ist. 


1. Aufstellung der Integralgleichung. 


Die Randbedingung am Loch erhält man für die Unterfunktion durch die Anwendung 
der Laplace- Transformation (8) auf die Randbedingung (7.1) zu: 


oU U ig 
en ET EZB ü a 2 ee en Fe Le 16). 
a Fo; für = a (16) 
Setzt man (15) und die Ableitung der Unterfunktion nach r ®) 
pr 
0 K,\— 
OU (r, p) (& ) E E p e) 
“3 N ae 22 24 rn 17 
dr A(p) dr A(p) „ Kı = +7 & c (17) 
in Gl. (16) ein, so folgt daraus der in der Lösung noch offen gebliebene Faktor 
Eh i 18 
ee B 
Pre R2 38 
und damit die Unterfunktion der Verschiebung 
le) 
uk B= 
BE BEE a RL TR ur 19 
DTRDETE 2) 9 x, Re 2 ra 


Bevor wir zur Übersetzung der Unterfunktion in den Oberbereich übergehen, wollen wir dimen- 
sionslose Größen einführen, um die Ergebnisse der Rechnung in möglichst allgemeiner, d.h. von 
der Randspannung q, der Abmessung des Kreisloches und den Materialkonstanten unabhängiger 


Form bringen zu können. Diese dimensionslosen Größen wollen wir mit überstrichenen Buch- 
Staben bezeichnen: 


REN RL a a er Bere 
a a C agq qa (20) 
9a ou our na 7 s 
ee en —, -=U =, 0,=—, = 
or erg’ 9 bog q q 


:)S. Watson, 8.202 und 203, Gln. (1) und (2) 
°) vgl. Watson, 8.79. 


5 


28 Re Air aus en) en kommt in den | dimensionslosen Größen | nur noch die Querkontrak- 
 tionszahl » explizit vor: TE EM ARE EEE TE er N 
# 5 az die Laplace- Melkstoemerte ER, a RE RER AR | 


5 ae “ war de N 
u RE VE IE | v | IE TE 
® erhält man aus ( “9 und ( @): ee ne ET Be RAN . 
Eine) I WET TE RENTE Be 1 SE “ er ee £ E zu EURE, 2 
© ea: eye N 
ee u6, = = ur DI 
Be e ErTeT ee ee 
Br: a Wir haben nun ae Ergebnis in den Überbh RT kranbe etzen Am nächsten 


0 läge es, die Oberfunktionen zu Z und 1/N zu bestimmen und dann’den Faltungssatz anzuwen- 
0. den ?).. Dann hätten wir die Lösung unseres Problems in expliziter Form. Es ist jedoch nicht 

gelungen, die Rücktransformation von !/N auszuführen, dagegen konnten wir die Oberfunk- er 
2 tionen von @Z und ee N finden. Man kann dann die Gl. (22) in die Form a 


x RR r 


N ee PREIS: Ul,D)erN( Pe @Zle.p). BE 
Er b 

E: darauf den Faltungssatz anwenden und so eine Integralgleichung zur Bestimmung von % ge 
= n 'winnen. Wir zeigen, wie das möglich ist. 2 f ; 
4 Zur Rücktransformation von e®N und erZ benutzen wir die S chläfl ische Integral- 
a 1 S | | 
Es darstellung für K, und = ? RER 

3 a 

£ - K,0, 9) = fe” —— dx 

“ - =: 7 2 — r2 - 
3 j 5 . (24), 


die für Re p >0 gilt. Setzen wir zunächst r=1 und @=t +1, so können wir N (p) sofort 
darstellen in der Form | 


x 


ne on 
) r 
1 
mit KA nfaHN?- age Te 
(en 
Damit ist die Oberfunktion zu erN(P) schon gefunden: Bi ; 
BD Na) lt) ee ee Ay 


Für die Rücktransformation von @L(, p) gehen wir aus von der Beziehung 
Er 1 EZ ae 2.172 
zen <- Kenn [er Ri, de = 
Re PROB, 
2 EI Er ee 
. — ax |. EEN 
li R = Ya —r? 5 In DT E [a 
die sich auf Grund von (24) leicht durch partielle Integration bestätigen läßt. Setzt man 
hier wieder <= t +1, so erhält man auf der rechten Seite der Gleichung ein Integral zwischen 


?) S. Doetsch, 8.161. 
s) 8. Watson, 8.172. 


den Grenzen von #—1 bis oo. Definiert man eine Funktion ER ., 
3 ’ 0 füro<t<sr—1 
wa NE m m EI EN ee 
ZEN Dr IE Hi? 5 In e > 
so kann man die Integration auch zwischen den Grenzen 0 und w erstrecken und erhält: 


Z(,$) = | Dar Dr REES (30). 


oder kurz . 
DR De N Eee (31). 


Damit ist auch diese Oberfunktion bekannt. \ 
Wir wenden jetzt auf die linke Seite von (23) den Faltungssatz an und erhalten als 
ihre Oberfunktion das Faltungsintegral: 


ri 
U (r,p)e®r N (p) u A Ele EEE (32), 

wobei nach Gleichung (26) 
ER en er. 
Ye +1 — 1 


ist. Die Anwendung des Faltungssatzes ist an die Voraussetzung gebunden, daß die Laplace- 
Transformation für alle beteiligten Funktionen absolut konvergiert. Das ist für die beiden 


Funktionen % (7, ) und & (t) bei Rep >0 in der Tat der Fall. Für die erstgenannte Funktion 
aus dem Grunde, weil sie für die Randbedingung , = —gbeir=abeschränkt?) sein 


muß, und für k (t), da diese Funktion bei t = 0 eine integrable Singularität hat (s. Gl. (26)), 
während sie bei hinreichend großen t > T wie (1— »)t anwächst. Der Anteil des Laplace- 
Integrals im Bereich T <t< © konvergiert aber für f(t) =? wegen Rep >0 absolut: 


a Fe u or 
BER —tRep Ber 
Rep ü ( 5 Re 5) | . 
de 73 
—- ze Br En 
Rep . T Rep 

. „Wegen der Eindeutigkeit der Abbildung durch die Laplace-Transformation müssen 
die Oberfunktionen der beiden Seiten der Gleichung (23) gleich sein, und aus (31) und (32) 
erhalten wir: 


. (33) 


Sy | 
[e>} 
| 
S 
= 
Er 
| 
IN 
—S 
IS 
\ 
e>) 
» 
Ss 
1 
Qu 
>] 
I 


[2 = = 
[Te Mkti—T)dır=z(r,t) re u rs (34). 
) 
Das ist in bezug auf % (7, ) eine Integralgleichung erster Art vom Volterratyp, die sich leicht 
numerisch auswerten läßt und die wir als endgültige Lösung unserer Aufgabe ansehen wollen. 


Der Kern k ((— 7) der Integralgleichung ist singulär, aber integrabel. 


2. Verschiebungen und Spannungen hinter der Wellenfront. 


Aus der Integralgleichung (34) können wir sofort eine wesentliche Erscheinung in der 


Ausbreitung des Stoßvorganges ablesen. Da für jedes € <?F—1 die Funktion z (7, d.=l ist 
(s. Gleichung (29)), so ist auch 


ar, de0 früheren (35), 
d.h., bis zut= 7 — 1, also bis ct = r — a erfahren die Scheibenpunkte r von dem Stoßvorgang 


noch nichts. Mit anderen Worten: die durch den Stoßvorgang hervorgerufenen Schwingungen 


haben eine Wellenfront ‚ die sich mit einer endlichen Geschwindigkeit c = YE’Jo aus- 


breitet (s. Gl.(5)). Derselben Erscheinung begegnet ma i Längs- EL 
gungen in Stäben }). BEBER n auch bei Längs- und Querschwin 


°9) 8.Doetsch, 8.14. ” 


en 
2 Tann ie Be a Are 
ee a ee 

a FE RS ze 5 BE =f( (, Dürr u Tre n 
nF = we Aus Stetigkeitsgründen muß die Funktion % (r, "2 ae nach dem Durchälkreiiegt 


. lee engen a 


| =ala- | Hemer Am + HEeF@l nat re) = 
Re seta et 


= 
Er 
2 
3 


der Wellenfront, d.h.fürti=r—1 re mit e— 0 gleich Null sein. Zur Kontrolle setzen wir E: 
in Gl. (86)t=r—1-+eund gehen dann zur Grenze e> 0 über. Da der Kern k (7) in dem 
eisiegrablatsbereich, das Vorzeichen nicht ‚ändert, so können „Wir nach dem NER He 
 Integralrechnung schreiben: es 


nt ag 


1lts 


ni 


"Nach Potenzen von & Er hackelt, erhalten wir tatsächlich En 
ER $ ne a en 5 = 2 2 OF E +-höhere Glieder FE 

UNE 7 Ur am en een re ER EN 52 

Er, u) V) 2. + höhere Glieder Be: 

- ar Berechnung der r Spannungen o,und 5, und der Teilchengeschwindigkeit @hinterund 3 

an der Wellenfront müssen wir aus der Integralgleichung für % noch solche für % und % ermitteln 
Durch Differentiation der Gleichung Bon nach r ergibt . bei Beachtung der Gleichung (37). a 
folgende ee für «: 
% 


+1 


17 (r,T) k d- 94-55 


1—1 


Ve+ arwan | 


HH | 


1 
+1 


Die Integralgleichung für % erhält man aus (36) durch partielle Integration. Führt man 
für das unbestimmte Integral des Kerns die Bezeichnung ein 


k*(t—7)= - [| Be er year tnazı 


ande EERhFLARE (39), 
zen [i-: Ze + i—T+ 11]! m e—: Sara Te =, 
so schreibt sich Gleichung (36) in der Form 


(v6 (r,T)dk*( t—N) = rd. 
tl 


Integriert man die linke Seite partiell und beachtet, daß % (r, t=r—1)= 0 und k* (0) = 0 
ist, so folgt daraus die Integralgleichung für %: 


N 2 CA Er (40). 

r—ı 
Sie hat dieselbe freie Funktion f, (r, t) wie die Gleichung für @ selbst, aber einen anderen Kern. 
Für die weitere Untersuchung ist es zweckmäßig, für % eine Integralgleichung mit dem gleichen 
Kern zu haben wie für @ und @. Zu diesem Zweck differenzieren wir die Gleichung (40) nach t. 
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Mit k* (0) = 0 und 


ok d— 3) ___ Okt ea rg NN EN AERRe, = 
ot RE 


erhalten wir dann die Integralgleichung für @ in der gesuchten Form 
; - 1 En _ = 
[we Hdrdddi=- Ze +Nior u Pa a 1 ee (41). 
r—1 


Auch in den Integralgleichungen (38) und (41) setzen wir für die obere Integrationsgrenze 
wie oben =? —1 + e und lassen dann e> 0 gehen. Dann erhalten wir wieder mit Hilfe des 
Mittelwertsatzes: 


+ + höhere: Glieder 
Te I a 
Y2 + höhere Glieder Yr 


y 2 y ; 

“(r,t=r7—-1-+0)=lim = + höhere Glieder 1 
| el]. 0.2000. (42.2). 
Y 2 e-+ höhere Glieder Yr 


Die ersten Ableitungen %’ und %' besitzen also an der Wellenfront eine Unstetigkeitsstelle. Setzt 
man (37) und (42.1) in das Elastizitätsgesetz (21) ein, so ergeben sich die Spannungen un- 
mittelbar hinter der Wellenfront zu: 


o,r,t=r—1+0) Be #-)# ET I (43.2) 

r r 
Für r=a und t= 0 erhält man aus (43.1) entsprechend der Randbedingung (7.1) o, = — 4. 
Wegen v=O0beir=aundt=0ist 0o,= vo,= — vg, was durch Gleichung (43.2) bestätigt 


wird. 

Da vor der Wellenfront die Scheibe spannungsfrei und in Ruhe ist, so geben die Größen 
°,, 5,und « für t=r— 1 +0 den Sprung in den Spannungen und in der Teilchengeschwin- 
digkeit, und die Gleichungen (43) und 42.2) geben die Gesetzmäßigkeit, mit der sich diese Sprung- 
größen ändern. Aus (42.2) und (43.1) erhält man noch das Verhältnis zwischen den Sprung- 
größen für die Geschwindigkeit und für die Spannung 9, zu: 


oder mit (20) und (5) 
aut 1 
mi w Vow na 45 


Bemerkenswert ist hierbei, daß dieses Verhältnis nur von der Größe des Wellenwiderstandes 
WE Yo E’ und nicht von der Art der eingeleiteten Störung beir = a abhängig ist, eine Erschei- 
nung, die auch von anderen elastischen Ausbreitungsvorgängen her bekannt ist. 

Wie man aus (42) und (43) sieht, nehmen die Sprunggrößen mit Ya/r oder allgemeiner 
ausgedrückt, mit der Quadratwurzel aus dem Verhältnis der Ringfläche 2x a ö (ö = Scheiben- 
dicke) an der Stelle, wo die Sprunggrößen eingeleitet wurden, zu der Ringfläche 2x r ö an der 
Wellenfront ab. Dieselbe Gesetzmäßigkeit hat Kucharski für Vorgänge an einer durch 
Längsstoß hervorgerufenen Sprungstelle in geraden Stäben mit veränderlichem Querschnitt 
gefunden !). Daß diese Übereinstimmung trotz der Ringspannung o,, die das Scheibenproblem 
vom Stabproblem im allgemeinen unterscheidet, stattfindet, ist darauf zurückzuführen, daß 
die Ringspannung zu den höchsten Ableitungen der Differentialgleichung (4) keinen Beitrag 


liefert. Diese höchsten Ableitungen bestimmen aber die charakteristischen Grundkurven, auf u 
denen die Unstetigkeiten existieren. 


wir die RER ” " oder 5 5, ‚oder 5 Gr 


7 1. as X ° r schen ‚Reihen. berechnen: 


Haöeet EN ET 6 1 +0 a 7 


„4 ”, „ii N”, 


ee zart +0) 4 6 1 +0 See 


TH nz Der +10, ie d+i-0, + 


- Da % 


(A Io (a An, 


| ea eu HH. 
Die erste a die ee für ein gegebenes 7.nach dem Durchschreiten der Wellenfront j 
> r —1), die zweite für ein gegebenes £ hinter der Wellenfront (r<t +1) ” 


"Wir stellen noch die Ergebnisse für die beiden ersten artiellen a it © 
er h er. 
von u‘, 2 und 6 5, an der er zusammen: en i unge 2 : = i ; 


: DR CE RE 
NE A ere NR RT hi) 
a a ei N 5 
FRE x | 647° Ar ie : | > 
„A r Fe 
i - mn, Be; 
Be: i eler st) ; 
£ Se a ee Er 
i 02 rlerars a 
e S ut an 237. Ha Fee 11» ale 7 11» a Ei 
; ; At ara euren) 
3 1 21. 9 7,» 7 1» 
a eh u Br a Dr BB ern 
s el. dar m ER PE ») 
7 ee 
u ' u Im ee a 2 
: RS ( Fre t”) © 
: er fl 7v 7» 2 
R- EL SH: ee ) 
| ; ee ee 49). 
a a n 1102 = 
; : gel r» 64r ee Ts Er: =) 
ash 27 0 Im, Hr, 
# = (er EEE ») 


‚Aus (42.2), (43), (47), (48) und (49) können wir die ersten drei Glieder der Reihen (45) und 
(46) berechnen. Damit läßt sich ein kleiner Bereich At oder Ar hinter der Wellenfront unter- 
E: suchen. Für größere Bereiche müßte man noch höhere Ableitungen an der Wellenfront her- 
A leiten. Da dieser Weg sehr mühsam ist, so verfolgen wir ihn nicht weiter. 
Für beliebig große Bereiche ermitteln wir die Zustandsgrößen der Scheibe durch nume- 
rische Auflösung der Integralgleichungen (36), (38) und (41). 


en a MEER 
j Eu 2 - or TORE u he EEE ar A | 
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_ - r ze u 2 ne Bor - ey Niger er 
Den Integrationsbereich von —1 bis t unterteilen wir in n gleiche Abschnitte & und ER 
führen die Bezeichnungen ein: | A en E: £ z 2 $ a, 
Der—Lherltht. eher ee 
H=r—l u=r—-1+6, mer—l +möd. En ; F re A - 


Damit schreibt sich die Integralgleichung für % folgendermaßen: | = 
er De E ER = IE 
I Sardeh- Harlem ehrn 2.2.2... 0 


a 4 } 
Nach dem Mittelwertsatz der Integralrechnung erhalten wir hieraus: 


R Zu | 3 
ED nn tm en mael= RR RN 


m=1 A 
mit 0 SNnm Sı. 
Die Integration unter dem Summenzeichen können wir ausführen: 
tm 
f kn —r)de= km 
im—ı 
19 14m m+1)) VOF mm + In [14m m+1) + 
RL 2 (52). 
SE He 1—r a 
+/U+% m 1) 9% —1] 5 (1 + (n—m) 8) / (l + (m —m) 1 — 
14» TE le RE 
— Zn [1 + m—m)ö+ (1 +R— m) —1] 
Setzt man das in (51) ein, so erhält‘ man folgendes Gleiehungssystem: 
ne x - 
aD (53). 
m=1 


n(n+1 
Dieses System von n Gleichungen ist noch exakt, kann aber nicht aufgelöst werden, da es | ai 


Unbekannte enthält (s. den Aufbau des vereinfachten Gleichungssystems weiter unten). Nun 
nehmen wir näherungsweise an, daß 72m von n unabhängig sei, und setzen 


Nnm — Nm | 
an BE, EP un nl ea, 54). 
ur, m—ı + Nm ö) — Um J 


Dann erhalten wir aus (53) ein System von n Gleichungen 


mit den » Unbekannten %„, welches angenäherte Werte für % (r, £) liefern kann. Die Näherung 
wird um so besser sein, je kleiner die Intervallteilung ist. Wir veranschaulichen noch den Aufbau 
des Gleichungssystems (55), indem wir einige Gleichungen anschreiben: 


Ju = % ku 3 
Jia = “ ka + u) kys 
Jıs = U kyı + Ug kb gg + U kzs 


Dieser Aufbau, der dadurch gekennzeichnet ist, daß oberhalb der Hauptdiagonale keine Glieder 
vorhanden sind, ist charakteristisch für Gleichungssysteme, die aus Volterraschen Integral- 
gleichungen entstehen. Da die Koeffizienten knm nur von der Differenz n—m abhängig sind 
(s. G1.(52)), so sind die Koeffizienten in jeder zur Haupt-Diagonale parallelen Schrägreihe unter- 
einander gleich, z. B. k,, = kg, = kn und %y, = ky, = ky, usw. Somit brauchen nur die in der 
ersten Spalte stehenden Koeffizienten berechnet zu werden. 


Bes 


or: 


Die Gleichungssysteme für @ und @ haben dieselben Koeffizienten knm wie das für %, 


® ‚da die zugehörigen Integralgleichungen (38) und (41) denselben Kern haben wie (36) und unter- 
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scheiden sich von Gl. (55) nur dadurch, daß dort f,„ oder f,, statt f,„ einzusetzen ist.' Außerdem 
sind die kn, von r unabhängig, die fin, fan und fzn müssen dagegen für jedes 7 neu ermittelt. 


. werden. Wir haben %, #, @ und daraus o, und ö, für r= 1; 2,5; 4; 8; 12 und 16 berechnet. 


Für das Zeitintervall ö haben wir im allgemeinen den Wert ö = 0,5 genommen. Zur Kontrolle 


haben wir für r —= 4 außerdem eine gröbere Intervallteilung gewählt (6 —= 1) und hinreichende 
"Übereinstimmung mit der genaueren Rechnung gefunden. 


Die Auflösung der Gleichungen (55) und solcher für % und % gibt uns die Größen 9, 6g: 


und u unmittelbar in Abhängigkeit von der Zeit t mit 7 als Parameter. Die Rechenergebnisse 
für » = 0,3 sind in den Bildern 2, 3 und 4 aufgetragen. Hierbei ist zu bemerken, daß die Größen 


Ir m usw. (m = 1,2,3...) irgendwo innerhalb des Zeitintervalls = — in ı = Öliegen müssen. 
Für größere m kann man sie in der Mitte dieses Zeitintervalls auftragen. Für kleine m dagegen 


(m = 1,2), d.h. gleich nach dem Durchschreiten der Wellenfront, kann das zu größeren Ab-- 


weichungen vom wirklichen Kurvenverlauf führen. Hier konnten wir aber diesen Verlauf in 
Zweifelsfällen mit Hilfe der ersten drei Glieder der Taylorschen Reihe (45) mit aus- 


‚reichender Näherung berechnen. 


Die Kurven in den Bildern 2, 3 und 4 haben alle denselben Charakter. Bis ut=r—1 
bleibt die Scheibe spannungslos und in Ruhe. Zur Zeit = r —.1 erreicht die Wellenfront, die 
sich mit der Geschwindigkeit c ausbreitet, den Kreisumfang vom Halbmesser 7, da der von der 
Wellenfront zurückzulegende Weg gleich r — 1 ist. Das Eintreffen der Wellenfront ist durch 


den Sprung, den die Spannungen 0, und 0, und die Teilchengeschwindigkeit @%° dabei erfahren, 


gekennzeichnet. Der Spannungs- und Geschwindigkeitssprung selbst läßt sich aus den Glei- 
chungen (42.2) und (43) errechnen. Nach dem Durchschreiten der Wellenfront ändern sich 
0,, 0, und %' ziemlich steil in entgegengesetztem Sinne zur Sprungrichtung und nähern sich 
dann nach einem Hin- und Herpendeln langsam dem stationären Zustand. Die Teilchengeschwin- 
digkeit %: strebt dabei dem Wert Nullzu. Der stationäre Wert von o,hat ein anderes Vorzeichen, 
als der an der Wellenfront, da zuerst wegen v=0 die Ringspannung 0, = vW =», = — vg" 


Se < 0 ist, während später infolge der inzwischen eingetretenen radialen Verschiebung « 


in der Umfangsrichtung Zugspannungen entstehen. 


“ Das Fortschreiten der Wellenfront ist besser aus der Auftragung über 7 mit t als Parameter- 


in den Bildern 5, 6 und 7 zu erkennen, die aus den vorhergehenden Bildern ermittelt wurden.- 


"Zu einer Zeit £ befindet sich die Wellenfront beir =t +1, da sie bei r = 1 ihren Ursprung hat. 


Auch aus den Bildern 5, 6 und 7 ist der Übergang zum stationären Zustand hinter der Wellen- 
front zu ersehen. 

Dieser stationäre Zustand der Scheibe ist im vorliegenden Falle dadurch gekennzeichnet, 
daß die Zustandsgrößen o,, 0, und “' konstante Werte annehmen (s. Bilder 2, 3 und 4). Dabei 
ist der asymptotische Wert der Teilchengeschwindigkeit %' gleich Null, d.h. die Scheibe geht 
in den Ruhezustand über. Dementsprechend muß die Spannungsverteilung in die für sta- 
tische Belastung der Scheibe durch den Druck g am Lochrand übergehen. Für diese Rand- 
bedingung sind die statischen ‚Verschiebungen und die Radial- und Umfangsspannungen, wie 
leicht nachzurechnen ist, : 


_ 1 
 (d—»)r 
1 
= ER (56) 
= 1 
SR Age: 


Die Kurven in den Bildern 2 und 3 nähern sich tatsächlich an die nach Gleichung (56) dort 
eingetragenen asymptotischen Geraden. Auch die Bilder 5 und 6 bestätigen, daß weit hinter 
der Wellenfront die Spannungsverteilung in die statische Spannungsverteilung übergeht. 


IV. Zusammenfassung. 

Wir haben die Ausbreitung von Wellen in Kreislochscheiben untersucht, die durch eine 
plötzlich aufgebrachte drehsymmetrische Störung am Lochrand ausgelöst werden. Dabei haben 
wir hier angenommen, daß auf den Lochrand plötzlich ein radial nach außen gerichteter, zeitlich 
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endlichen Geschwindigkeit c= 


er zu dem Spannungssprung von dem Wellenwiderstand W = 7 
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| konstanter Druck q wirkt. Es hat sich geze st 
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das Verhältnis des Geschwindigl 
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Wellenfront geht die Scheibe in einen stationären Zustand, im vorliegenden Fall in den R e 
zustand, über. Dementsprechend geht die Verteilung der Verschiebungen und Spannungen in 


- 


diejenigen Spannungen und Verschiebungen einstellen, die eine statisch erzwungene Randver- 


schiebung der augenblicklich vorhandenen Größe u = v, t hervorrufen würde. a 


Eingegangen Oktober 1947. 
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Über eine Verallgemeinerung des Verfahrens der Kombination 
von Newton’scher Methode und Regula falsi zur Auflösung 
| einer Gleichung f (2)=0. a 


Von Erwin Pflanz in Marbach/Neckar. Bu 


In Verallgemeinerung des kombinierten Verfahrens Newton-Regula falsi wird gezeigt, daß sich eine 
(reelle) Wurzel x&* einer Gl. f(x) =0 durch die Nullstellen gewisser Paare von Polynomen N-ter Ord- 
nung (N >2) Hermite’scher Interpolation einschließen läßt. Hierbei wird insbesondere Monotonie 

_ der N-ten Ableitung von f (x) im betreffenden Interpolationsintervall vorausgesetzt. Bei wiederholter Durch- 
führung konvergiert das Verfahren gegen x*. , 

The combination of-the method of Newton and the „regula i* is generalized, and it is demon- 
strated that a (real) root x* of an equation f (x) = 0 can be included by the zeroes of certain pairs of (Her - 
mi te’s) interpolation polynomials of the N-th order (N > 2). Itis especially assumed that the N-th,deri- 
vative of f (x) is monotonous in the interval of interpolation concerned. By iterating the process the method 
converges against &*. 

Ilpn o6oömennn MeTona KOMÖHHHPOBAHHOTO IIPHMEeHeHHA CIOCoßa HEWToHa H CTOCOoda 
IPOCTOrO HHTEPINOANPOBaHun (regula falsi) aBTop IOKASBIBAET, YTO (BEINECTBEHHEIÄ) KOpeHR x* 
ypasHeHna f(®)=0 MoMeT ÖBITB 3ARIIOyeH Meikıy KOPHAMH HEKOTOPHX Map HHTEPNOAIA- 

. HMOHHEIX IOJIUHOMOB IpMuTa N-Horo nopanka (NZ 2). Ilpu oTom B yacTHocTu npen- 
NONOTAEeTCH, YTO N-Haa IPOonsBonHaA OT [(%) B PaccMaTpmBaemoi OÖNACTH HHTEPNOAANHH 
IOCTOAHHO BO3PACTaeT MIIH IHOCTOAHHO yöslBaer, Ilpm NOBTOPHOM NPHMeHeHHH 9TOTO 
mpomecca, 9TOT MeTON IPNUBOAHT K pemenmo ®*, 


1. Einleitung. Um die Wurzel «* einer Gl. f(x) = (0 numerisch festzustellen, kann 
man etwa davon ausgehen, f(x) durch eine Funktion p (x) zu ersetzen, welche an m Stellen 


2 u 
# 


.. 


z&,(u=1,..,.m) in der Umgebung von &* mit f(x) in den Funktionswerten und den aufein- _ 


anderfolgenden Ableitungen nach & bis zur (a. — 1)ten Ordnung (a, —=2,3...) übereinstimmt 
und dann als Näherung für x* eine Nullstelle von 2 (x) zu nehmen. Die Ersatzkurve y= p (x) 
berührt also y= f(x) in den Punkten x = x, in (a, — 1)ter Ordnung, ihr „Annäherungsgrad“, 
d.h. die Gesamtzahl ihrer mit f(x) gemeinsamen Funktionselemente (einschließlich der ge- 
meinsamen Ableitungen) istr=a, +... +0%,>m. Praktische Vorteile wird. diese Methode 
nur dann aufweisen, wenn die Nullstellen von p (&) leichter zu berechnen sind als die von f (x) 


und wenn der Fehler der damit für ©*erhaltenen Näherungsich nicht zu schwierig abschätzen ° 


läßt. Man erhält so bei Annäherung durch Polynome (nr — 1)ter Ordnung P? (x) speziell für 
n = 2 bei zwei Interpolationsstellen x, x, die Regula falsi, bei einer einzigen Interpolations- 


“ stelle x, das Ne wto.n’sche Verfahren. Für m = 1, n >2 sei insbesondere auf die erweiterten 


Ne wton’schen Formeln hingewiesen!) 


!) Näheres über die erweiterten Ne wton’schen Formeln siehe etwa: Fr. A. Willers Methoden der prak 
! VERA, Ss 
Analysis, 5.175, Berlin 1928; G. Schulz, Formelsammlung zur praktischen Mathematik, 8. 48, ne a 
Berlin 1937; E.Bode wig, Comment. math. helv. 8 (1935/36) 8. 1—4. Die Annäherung einer Funktion in der Umgebung 


een durch einen Kegelschnitt in allgemeiner Lage untersucht R. Beyer, Z. angew. Math. Mech. 16 (1936), 


u ee a u a 


(siehe Bild 1). Mit den Nullstellen gewisser 
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Beibehaltung der Zahlwerte der %, mit «12, x12] 
als neuem Ausgangsintervall usw.) entweder 


.aber die Folgen der Näherungen x], x] für 
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Yale ıH NISEnden nn nun der Fall m = 2, n—= 3, £* an Hierbei den wir 
in Abschnitt 2 (bei festem n) Paare von Polynomen pe), ? (x) vom (höchstens) (n — 1)ten 


Grade angeben, welche unter Voraussetzung der ot der (a —I)ten Ableitung fa, ” 


( 
im betreffenden x-Interwall f(x) ‚einschließen 3 


dieser Einschließungspaare.erhält man dann 
Näherungen «2, x!2] für x*, die zugleich «* _ 
einschließen und die besser sind als die Aus- 
gangsnäherungen x, = all, = all. Dieses 
Verfahren ist ersichtlich eine Verallgemei- 
nerung der bekannten Kombination 2) von 
Newton’schem Verfahren und Regula falsi 
(siehe Bild 2). In Abschnitt 3 zeigen wir, daß 
bei seiner wiederholten Durchführung (unter 


unser Verfahren nach endlich vielen Schritten 
exakt x*liefert, also von selbst abbricht, oder 


v— a von beiden Seiten monoton gegen «* Bild 1. 

streben. Dabei ist vorausgesetzt, daß im Innern 

des Intervalles (x, x;) f(x) nur eine einzige Nullstelle besitzt. Schließlich werden in Teil 4 für 
n = 3, 5, 7 Formeln zur Gewinnung von Näherungen für x* auf obigem Wege angegeben und 


_ diese auf Beispiele angewendet. — Gegenüber dem Ne wton’schen Verfahren und der Regula 


falsi bzw. deren Kombination hat man mit der vorstehend (für n>3) skizzierten Methode 
beim einzelnen Schritt eine größere Rechenarbeit zu bewältigen. Insbesondere sind hierbei 
für n > 3 Ableitungen zweiter und höherer Ordnung (deren Existenz und Stetigkeit voraus- 
gesetzt) von f(x) an den Intervallenden zu berechnen. Andererseits kann man so bei gleicher 
Intervallänge 2= x,— x, > 0 mit einem einzigen Schritte x* näher kommen als im 
Falle n = 2. Zugleich ergeben sich dabei obere und untere Schranken für x*, so daß sich eine 
besondere Fehlerabschätzung erübrigt. 


2. Aufstellung von Einschließungspolynomen und Erläuterung des Verfahrens der 
Wurzelbestimmung. a) Die Funktion f (x) und ihre «-ten Ableitungen f(®%(x) nach & sollen in 
den m-Interpolationsstellen &, (u= 1,2, Mm) die Werte 


| Ta) =f0 a), Fa)... fd (2) 
besitzen. Es sei 
: m > 1, 1, 2 0m 1: 1 Hd tr. Fon N. 


f(x) seiin dem betrachteten x-Intervall eindeutig und bis zur r-ten Ordnung stetig nach 
x differenzierbar. 

Es gibt bekanntlich ?) genau ein Polynom y= Pa, ...a, (2) vom Grade N<n—1, 
das in den Punkten &, mit f (x) sowohl hinsichtlich der Funktionswerte als auch der Ableitungen 
nach x bis einschließlich der Ordnung &, — 1 übereinstimmt. Wir denken uns für das folgende 
dieses Polynom bestimmt. Dann gilt ®) 


: 1 PR 

| fi) = Pa,, 0 (%) + rn (%) = Pas, ... m (%) A m (29) ee zu) We} 

In (1) ist z, eine gewisse uns hier näher nicht interessierende Zwischenstelle im Innern 

de Intervalls A<z< Bmit A = Min (z.. ... &, 2), B= Max (z,, ... 2%, &). Aus (1) ergibt 
sich speziell bei zwei Interpolationsstellen &,, % 


1 
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2) Vgl. Fr. A. Willers(a.a. 0.1), S.171 und A. Loewy, Pascal’s Repertorium der Höheren Mathematik I, (1910), 
S. 354, sowie die dort angegebenen Literaturhinweise. 

s) Vgl.etwa A.Markoff, Differenzenrechnung (deutsche Übersetzung von Th. FriesendorffundE.Prümm), 
S.1—3 und 8.6—8, Leipzig 1896. 

4) Eine andere Darstellung des Fehlergliedes r, (x) in (1) und zwar durch einen geschlossenen Integralausdruck findet 
sich bei G. Kowalewski, Interpolation und genäherte Quadratur, 8. 37, Berlin 1932. 


FI) = Pa,, a, (%) a 


g* 


Y 


E ersichtlich auch aus (2), wenn man in (2) formal &, 


Übereinstimmung im Funktionswert und gewizßen aufeinanderlolgend ei 
‚gesetzt ist. 


ei gilt. Eine der beiden Ungleichungen (4a), (4b) trifft nach (2) sicher zu, wenn in =” <a = Bz 


Fin a) bis Ko En RER die Zwise nstell 


bzw. ?0,«, (x) die Polynome 9%, (x) bzw. Me rsteht. 
zugrunde legen werden, für 0,6% 20, u +8 =n2 t, 
den Fall, daß nur im einen der beiden Intervallendpunkte ER %, für f(x) u 


. y 
Unser erstes Ziel ist nun, bei &, og >0,% Pe ih utg=u> Pr En 4 A = 
Paare von Polynomen ?%,,a,(%), Pa,,5,(®) zu finden, welche in Eye Intervall Le z< u 
einschließen, so daß also dort stets entweder e: ei: .e- 
Paun (2) =! GES EICH ” ” ” “ EEE - > EUER 
oder E | ae SE be 


Dt 


ee ee . Bus - u EN. ehe 


(— 1): fiat) (z,) AST ER IR Er e 5) 


ist für die zwei oben erklärten Zwischenstellen 29, 2, (%ı < 2, 2% < %2)- Beschränken wir uns 
auf sts dr +%=n 22, so folgt aus (9) mit &, — 0 = — (&, — )=—4 N 
ae If 0) (x) na La FE RA ng 5 N Re VE 5 


(6) ist erfüllt, wenn A= +1, +3,... ist, zim Intervall z,< x< x, liegt = dort stets 
entweder | Be Br. 
2: | ae DT DEREN, Fe OR 


oder Ber“ 


Er er eu 2 


ist, d.h. wenn fd (x) in 2, <x<x, monoton ist, wobei f(” (=) dort stetig vorausgesetzt sei. a 
b) Wir werden nunmehr solche mit (6) sich ergebenden Paare p,,, (2), Px,5,(2) vn „Ein 0 

schließungsparabeln“ für n=2,...7 angeben. Dabei sei jeweils eine der Ungl. (7a), (7b) als Bir 

erfüllt vorausgesetzt. Durchmustert man sämtliche nach (6) in Frage kommenden Möglich- 

keiten (%&,, &), (&%, &%), so erhält man die folgende Zusammenstellung: 


Auf die Angabe weiterer Paare (bei n>7) sei verzichtet, da die Gesetzmäßigkeit bei 
deren Bildung erkennbar ist. Von Bedeutung für uns sind hier, wie aus 2c) hervorgehen wird, 
Einschließungspaare Pa,,a, (%), P&,,a, (®) mit &, & > 0,%, 5> 0. Ersichtlich gibt es für jedes 
n>3 mindestens ein solches Einschließungspaar. 3 
Der Falln=2, der auf die in der Einleitung erwähnte Kombination von Regula falsi 


und Ne wto.n’schem Verfahren führt, ist durch Abb. 2, der Fall n = 3 durch die Bilder 3-6 
veranschaulicht. 


c) Es werde jetzt für m = 2,n>3 das Verfahren der Wurzelbestimmung unter Verwen- 
dung der Einschließungsparabeln erläutert. Wir setzen dabei voraus: 


(A,) Es sei f (a1) f(m)< 0. 
(42) f(x) besitze innerhalb J (x, < x < x,) nur eine einzige Nullstelle «*. 


(B}) Die Polynome 2,,.,() bzw. p%,3,(x) berühren /(x) in 2, a,-bzw. %,-punktig, 
in %, &-bzw.&,-punktig. 


Ju u= — Ber —&g) sei eleich einer der Zahlen‘ ET, 1, ah, 


© Linz e< %, sei fm (2) stetig und entweder ständig < 0 oder = ne je-D er 
dort monoton. 


= Fa ” 

Wegen (B.), (B,) lassen sich Para (x) und 5 IE als Pol u 
N Zr ynome höchstens Be ten 

"Grades in x eindeutig festlegen. (B,) bis (B, ) und (6) zusammen gewährleisten nach (6). na 

u B) die u leßung, also das Zutreffen einer der ern (4a), AD) Re: = 
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Bild 4. Bild 5. Bild °6. 


-Pa,0,(%) und p%,z,(%) besitzen wegen P«,a,(%1) ' Px,,a,(%g) = P%,,3, (%1) ' Pa,,%, (2) < 0 
Nullstellen- innerhalb J. Für diese Nullstellen 2=1,(e=1,2..r) von 94,a,(2) und 
z=%(0o=1,2,...8) von 9%,x,(%) gilt, wie man leicht einsieht, entweder 
\ TE a RT EN = (08) 
oder 

ee a ET a ER ra > (8b). 


Man hat so mit t,, T, für x* Schranken gelunden; die jedenfalls näher bei &* liegen als x, 

und &,. Ist insbesondere für eine der Stellen u =% bzw... =ts ff) = 0 bzw. U) S=nlE 

- soist damit auch x* gefunden ; denn wegen der Voraussetzung (A,) ist dann Te = ar bzwe Hi ee, 

Bei den in 2b) angegebenen Polynompaaren, für welche nicht alle Zahlen &}, &, &1, &%, von 

Null verschieden sind, könnte unser Verfahren versagen, da dann u. U. die Einschließungspoly- 

nome innerhalb J keine Nullstellen besitzen (siehe etwa Bild 3). Mit unserer Annahme 

(B,) &, &, &, &> 1haben wirletztere Möglichkeit jedoch ausgeschlossen. Die für unseren Zweck 

geeigneten Paare (%,, &), (&, %) sind in der Zusammenstellung in 2b) durch Sterne besonders 
hervorgehoben. 


3. Konvergenz bei wiederholter Durchführung des Verfahrens. Zum Konvergenznach- 
weis wollen wir dasjenige der beiden Polynome, das in J unterhalb f(x) verläuft, mit ?,,.,(%) 
bezeichnen. Es gelte im folgenden also die Ungleichung (4a). Weiter darf ohne Einschränkung 


der Allgemeingültigkeit 
Be el la). (9) 
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vorausgesetzt werden. Dann gilt von den Ungl. (8a), (8b) die letztere. f un n % 
Eine Wiederholung des in 2c) angegebenen Verfahrens kommt nur dann in Frage, nn 25 

in (8b) die Gleichheitszeichen nicht zutreffen. en | 
Es sei nun z[2 eine der Nullstellen £,, «1? eine der z,. Dann ist, nach (9), 
Lad) <0, Le >0,- rare BT 
Ja (al <a< al) sei als neues Ausgangsintervall gewählt. Behält man us Werte 
&, Xp, üy, gdesersten Schrittes bei, so gilt für die Einschließungspolynome Pe, (k), p3.5.(2) des 
zweiten Schrittes in J[2] (wie mit (2), (4a) und Vor. (C) in Abschn. 2c) folgt) ar 


DE (a) < fe) < paa), Wade ar 
Wegen (10), (11) und Vor. (A,) bis (0) in Abschn. 2c) sind somit auch für den zweiten . 


Schritt die Voraussetzungen erfüllt, die den Überlegungen beim ersten Schritte zugrundegelegt 
sind. . R 


U n n 
FEBEIFR ENT. 


p,(«) bzw. pr (®) besitzen also Nullstellen 1” bzw. 2” im Innern von J!2; es 
sei =[®) eine der z[2}, «3 eine der el2l. Dann gilt nach (&b) 
<a <eEZ «<a <, 
und entsprechend schließlich nach (v— 1) Schritten für entsprechend gefundene Schranken 
«ie, ze] von =* . 
nz dä. SEN HIT A ee (12). 
Wenn für ein v in (12) zum ersten Male ein Gleichheitszeichen gilt, ist x* nach endlich 
vielen Schritten gefunden. (Dies tritt beispielsweise dann ein, wenn f(x) mit einem der Ein- 
schließungspolynome identisch ist). Falls in (12) für jedes noch so große » nie das Gleich- 


heitszeichen sich ergibt, streben die Folgen x], z["] gegen «* für v—o, wie wir jetzt zeigen 
wollen. 


Die Nullstellenfolgen x[®?} bzw. x], die zu den Polynomfolgen Pe sie) bzw. pl", (x) ge- 
hören, genügen dann nach (12) der Ungleichung 


m<all<...<al<ar<al<..<lcn EN er (13). 


Da beide monoton und be- 
schränkt sind, konvergieren sie gegen 
Grenzwerte, die wir mit ==] bzw. 
x{®] bezeichnen wollen. (13) ergibt 


Aaeri<r<zii) . ,d(. 
Wir werden zeigen, daß in (14) 
nur die Gleichheitszeichen stehen 
können, daß also 
er a a a ER 5 
ist. 

Wir bemerken zuerst,®'daß das 
Polynom Fe. (2) (welches ja in 
PU die Kurve y= f(x) a,-punktig, 
und in =" «,-punktig berührt, also, 
wie wir sagen wollen, dem Intervall 


PRd 7. („P=D, zp-N) „zugeordnet“ ist), für 
jedes v= 2,3... mindestens eine Nullstelle x!®) besitzt, für welche 
ale1 <eli< a ee ee ee er N: (16) 
ist. 
Zum Beweis von (15) nehmen wir an, es gelte 
al®l Sat RT Ne (17). 


Es sei p[®] (©) das dem Intervall («[*], x{®)) „zugeordnete“ Polynom. Dann gibt es ein 
A,>0,so daß für alle Punkte x des Intervalles I <x<al®] +4, 


5 1 
> ee) SO (18) 


ist (f(a{°]) ist positiv wegen «[°] > x* nach der Annahme (17)). 


an 


REN 3 AUS ud; 


ka 55 


ee, ae +0 


u 
TE Aus Stetigkeitsgründen läßt Sn eine Zahl 3 2 0 so angeben, daß Be alle 0 < ö,, Le 5 er = 
Fe und alle » in un ı< ae die Vene, ET ee 2% \ 


: ee . 2 2 = x ae Y ee (a) = Bee) < 1 pie). ar Pe 2 le a9) k 


ende. Bi al ar uns nt ia. 1” 
Y PATEE 


Dr ee 


bezeichne weiter ersaay. dasjenige Poinan a Ar in dei ä;die Kurve 
in «l®1 +6, a, Be berührt (dr 6,>0), das also dem „Iutergalls 5 


Yy = (x) %-punktig, 
2) zugeordnet Ist k- 


.f 


S 


u, £ 


st: Wir Hehken uns A, und ö ad den Ungleichungen (18), (19) bestimmt und im Fe: 


genden festgehalten. 


Au (18), 19) erhält man insbesondere \ er Far 
Be aile|oı,m|- I ums ra|>0 Ne 


für 0 = ö,, ,<6 und alle x des Intervalles alrl<asal®r] +4, 


‘Wir betrachten nun die Folge der Polynome pr 2 (@)- Wegen al > E38 und a! re 


für v—o wird für ein genügend großes v — v a1 j in das Intervall (2]*] — ö, x1=1) und «1 
sowohl in das Intervall (2{”], [1 +8) als auch in’das Intervall (a1=], al=] + A,) fallen. “Für 


das Polynom el (x) , das zur Menge der Boa 171»? (x) gehört, gilt dann u Arie 


| peu (x) )| = 0 “für az t< a, 
Dies ep reh aber der in Gl. (16) zum Ausdruck gebrachten aber daßjedes 
Polynom ES (x) zwischen x[*] und &{?-1] mindestens eine Nullstelle besitzt. Die Annahme 


(17) führt also auf einen Widerspruch, desgleichen (wie man genau so zeigt) die Annahme 
ar! <a*<al®l, Damit ist (wegen (14)) (15) bewiesen. 
Ersichtlich läßt sich vorstehender Beweis auch auf den Fall n = 2, also auf das in der 


‘- Einleitung erwähnte EumDEe Verfahren Ne wton- negula fals:! übertragen, wenn man 


voraussetzt: 
(4°) Es sei f (a1) (a) < 0; 
(B’) x* sei die einzige Nullstelle von f(a) in, <xz=<%; 
(C’) f(x) sei in 2, <x< x, entweder durchweg positiv oder durchweg negativ; 
(D’) Ausgangspunkt der Newton’schen Tangentennäherungen sei derjenige Inter- 
vallendpunkt x, (u= 1 oder =2), für den f (z,) *f” (=,) > 0 ist. 


Die (zusätzliche) Voraussetzung (D’) gewährleistet, daß die Endpunktstangente die x- 
Achse im Innern des Intervalles schneidet. Letzteres wäre sonst, wegen & - &°%,°%, = 0), hier 
nicht gesichert (siehe hierzu Bild 2). 


‘4. Formeln zur Gewinnung von oberen und unteren Schranken für x* bein>3; 


Beispiele. a) Die Aufstellung der Gleichungen der Polynome y= ?,,,.,(%) bietet im Einzel- 
fall keine Schwierigkeiten®). Wir geben (im Hinblick auf unsere Anwendungsbeispiele) diese 
Gleichungen nachstehend für ,=2, %&=1 sowie u„=1, 98 =2bein=3und fürs, =3, 
9 = 2 sowie = 2,%=3bein=5an. Zur Vereinfachung führen wir hierbei die Variable 


1 R ; 
=t1—— 5 (21 +2,) = 2 — 2 ein. Y= Pa,o,(%) bzw. y=f(x) werden wir in Abhängig- 


keit von 2 mit P,,..(©) bzw. F (©) bezeichnen, die x* entsprechende Nullstelle von y=F(£) 
mit 6* = 2* — am: 


Es bedeute 
=Fh)=f(@), /-=FemM=fa) 2 =f9 a), I=I9 a), 6=1,2..). 
Weiter si2=,— 2, >0. 
Bei der Darstellungsform - 
£ ” n—1 
& 1, Os E Y 
Pu) = Z(z) 


5) Eine explizite Darstellung der Gleichungen der Polynome für m = 2 und beliebiges n = a, + ©, ist bei A. Mar- 


- koff,a.a. 0.3), $S. 6 angegeben, für allgemeines m und n (wie Verf. aus einer Besprechung im Zentralbl. f. Math. 3 (1932), 


S. 267 "entnimmt) beiP. Johansen, Skand. Aktuarietidskr. 14 (1931), S. 231—237. Formeln für ?,,(2) bzw. 9.,3(#) 
inNewton’scher bzw. Lagrange’scher Darstellungsweise finden sich bei G. Kowalew ski, a.a. 0.4) S. 19—20 


bzw. 8. 29—31. 


Kombinatio en 


ergibt sich für die Koeffizienten cv) peispielsweise: En 
Auf, +3f.+2hfL 0 dm +3fr —2hfl 


2 cal = fr —f_ ct) = e2,1) ERICH r (21) = x Be 

4 —f, —f_—2hfl 49 = f_—f++2hflı x 2 

1680 —5f, +llfL—2Hlf; —AfL) +2MfL 169 =5f_ +11f; +2h(f2—A+) H2RFE) ER 
4 =) —h(f4 4) 29 — 00:9 Be i 
Elf —f_)—6hfl—2hfZ Bel —fr) +6hf —2Mff ü 
4 (fr) +hlf4 +) | 2 


1689 = —3 (4 —)+2h(f4 +2) +2RefL 16089 = —3(f-— fr) —2h (fl +24) +2FE 
Die Schranken für x*bzw.C* erhält man durch Bestimmung der in das Intervall u<e<a 

bzw. —h<t&<h fallenden Nullstellen der jeweiligen Paare von ‘Einschließungspolynomen. 
Man zeigt leicht, daß bi n= 4 +%>4 %'%>0 und f, = f’ (z=*) +0 für genügend 
kleine h P4,,.,(2) genau eine Nullstelle in x, < x < x, besitzt, während Letzteres bei n = 3, 


&°0,> 0, also bei 9,5 (x) bzw. P,,ı (X) wegen ?,,. (2%) * Pu (8) < 0 bzw. P5,1 (%1) * P2,1()< 0 
offensichtlich für jedes A der Fall ist. 


Diese (einzigen) Lösungen von 
P,,0) 0: bzw: Du) DE (23) 
bezeichnen wir mit 5 . 
Gr &: bzw. 2,0, = Im +22.0 . 
Dann gilt für die Abweichung der Schranken für x* (f’ (z*) + 0) bei kleinem h 


h3 714 % 2 
ia= tat 2) | +... ee, (24a) 
* 
und bei hinreichend kleinem h | e 

j h5 & \2]2 
GE Er Bess ER SE 

t:32 | 

Y A) c#.\2]3 
u niun 2ng Bess en 


Die weggelassenen Glieder sind jeweils von höherer als n-ter Ordnung in AR. In (24a),... 


* 
bedeutet f? = f@ (x*), weiter ist dort E | <1. Die Ordnung in A der Abweichungen 


h R 
21 — Ci... ändert sich nicht, wenn man auf der rechten Seite von (24a),... die Quo- 
tienten f'/f’, f/f’,... anstatt für «* an einer anderen, im Intervallinnern gelegenen Stelle, 


etwa #—= x, nimmt und dabei f’(x)+0 in x, <x < x, vorausgesetzt werden darf. 

Die Formeln (24a),... sind nützlich als Rechenkontrolle oder zur Vorhersage der Größen- 
ordnung der im jeweiligen Beispiel zu erwartenden Abweichung. Man erhält (24a),...unter Heran- 
ziehung der Bestimmungsgleichungen (23) für £$,,„, bei Reihenentwicklung der darin auftreten- 
den Koeffizienten. 

b) Die numerische Auflösung der Gleichungen (nr — 1)ten Grades (23) erfolgt, falls eine 
Rechenmaschine zur Verfügung steht, am einfachsten wohl durch Iteration und zwar auch für 


n = 3. In letzterem Falle ergibt übrigens die übliche Auflösung der quadratischen Gl. nach 
Umformung ®) 


* 2/- 1 I Jr 275 ; 
FE Io A h— Be ae 5" a 1= Sc a ee R 
& Feng Ns ar = Er0)am 
* 2J+ 1 Jam 2r ‚ 
C 3 h aa Vz a —— = _ n EIER ° R % 
12 1 ae (&ı2 GA; Ir f, + 0)(25%) 


Es gilt vor der (positiv zu nehmenden) Wurzel das + bzw. — Zeichen 
in (25a): je nachdem 
| (+ —- Fa >0 Dane 
°) Den Gl. (25a), (25b) ähnliche Ausdrücke für N äherungen 2 „x,von «*erhält man bei Heranziehung der Polynome 
Pu Re. 


ji ji . Auf diese Formel ist hingewiesen bei Fr. A. Willers, 
+ 
e), 


* 
P9,3 (X), 93,9 (9), z.B. 20,3 = u, — 


+ 
Z. angew. Math. Mech. 18 (1938), S. 200, Fußnote 


An hr >0 -baw. = ° 


Schließlich ‚sei Ense (ohne Herleitung) a daß im Falle n— 3, falls 


VERAU REIT EREN AIR a A 5 a A er a u 
u . \ Y\ n R 


vom trivialen Falle c(D — c(1,2 — (0) genau einen 


er Ba (fahr +2hfl) a = 2 a 
tt für al = ar e neben der Entwicklung (24a) die Ungleichung er 
RER nf +) © 4 
T a leı—atelsh- I. se FEAT ET TER 27) 
= Man 
gilt. a 7) wie hernach in 1 (28 a), (28b) it das Gleichheitszeichen nur dan (2, Ne Er 
‚und c(1:2) heide gleichzeitig null sind. eh = FB = = 
x Weiter ist bei eaniekeit von (26) | a . % Di: 
n h®_ M, Be. 
| 21 — ee a RR a ee a > (28 a) BL - 
Er | | 3 
(ter-E +0 Mi— Max |” (ojlina,= Ben] nt Be 
und bei f(«) +0inz<a<e, % | | B 
we 0 
21 el ae N CE 
(m, = Min |’ (2). Oma <tem,). ae ER a 
"Die Vor. (26) ist erfüllt bei Kurven’ y = f (x), bei welchen die von den Parallelen 2 2 | 


und 2<=x, ausgeschnittenen Abschnitte der Kurventangenten in den Intervallenden auf 
verschiedenen Seiten der Kurvensehne liegen 2 
(s. Bild 8). Wegen der bei Anwendung von (27) | 
bis (28b) vorausgesetzten Monotonie von f” (x) 
besitzen diese Kurven in 2, < x < x,, abgesehen 


Wendepunkt, wie man sich leicht überlegt. 


e) Wir wollen nun noch zwei ee 
beispiele geben. 


Die Funktion f (x) = e®— 3 «? besitzt genau 
eine Nullstelle x* zwischen #2, = —0,5 und 
2%, = —0,4. f’(x)ist monoton n— 0,5 Se <—0,4, 
die Parabeln P3, (&) und Pı. (£) schließen dort 
also f(x) ein. 

‘Nach Gl. (24a) ist die Abweichung 


. #1 — #1, = —1,6.10% zu erwarten, wenn wir 
‚bei f, f” anstatt «* den Mittelwert 2, = — 0,45 = z 
nehmen. Be 
Man erhält nach Berechnung von f_,f/,f+,f4 aus (2) die Gleichungen 
Bi _ 0,180 600 54 . .:+0,04024046 .. eo 
1, x CF, 2 
= —0,18044131... .+0,040 08123... en ; 
Deren Auflösung durch Iteration führt dann auf | 
Ä 51 = — 0,458 965...< 2*< — 0,458 957... .— 21.2. 
Mit dem gewöhnlichen Newto n’schen Verfahren erhält man ?), wenn man von der 
Näherung x! = — 0,5 für x* ausgeht, nach 2 Schritten «[3) — — 0,45 897, wobei eine Fehler- 


abschätzung | «18! — x*| < 7,8.10% ergibt. 
Das Beispiel BE ()=e — 2, bei dem sämtliche wen in jedem endlichen Intervall 
. mit wachsendem :n. 


Mit dem (verhältnismäßig großen) Ausgangsintervall a 0, u Lund J-=--1LIh=e-—2 


")G.Schulz, a.a. 0.1), 8.47. 


192 


a Der a ei 


Kleine Mitteilungen | 


erhält man als Annäherung für die in diesem Intervall liegende Wurzel &* = In 2= 0,693 14723 


nach einem Schritte bei Heranziehung von 
P3,1(@), Pıe(d,n=3: 
P3, (8), Pa2(d), n = 4; 


P32 (0), Pas, n=5:- 


0,673 < x* < 0,703 
0,6905 < x* < 0,6948 
0,6930 < x* < 0,6933. 


i ist di i 1 310% 
Mit P,s(l) und P3, (€) ist die nach (24c) zu erwartende Abweichung |24,:— 23. = 3- h 
Dasch. en von Regula falsi und Newton’schem Verfahren (n 2 ergibt 

a , Fan 


sich hier, wie zum Vergleich angegeben werde, 


1 
e—|1 


=0,88...<ar<0,73...=- 


= 


Hinsichtlich Abschätzungen über die Abweichungen der Schranken von z* im Falle n = 2 


sei auf die Literatur verwiesen ®). 


8) Siehe etwa: Fr. A. Willers,a.a. 0.1), 8.172; G.Schulz,a.a.0O. 1), 8.46; A. Loe w y, Pascal’s Repertorium 
der Höheren Mathematik, a.a. 0.2), 8. 352—354; J. Lüroth, Vorlesungen über numerisches Rechnen, S. 175—177, 


S. 193—194, Leipzig 1900. 
Eingegangen am 15. November 47. 


KLEINE MITTEILUNGEN 


Praktische Interpolation höherer Ordnung. 
1. Direkte Interpolation. 


Bei Vernachlässigung vierter und höhererDifferenzen 
von & ist!) 


Ö2e” = 2 en 62 ei? Tl w? et. 
Damit ergibt sich für die Funktionen 
e”,;, e”, coshx, sinhx 
die einfachere quadratische Interpolationsformel ?) 
u 
| + %)=(1—u) (1-0 2) %tuyıl (1,1) 


und die einfachere Everettformel für kubische Inter- 
polation °) 


I 
S®+h) = s|sw-ua4s) @°] Yo 


E (1,2). 
+u[6—s (1+ u) @®] 9, 


Die entsprechenden Formeln für cos x und sin x gehen 
aus”diesen hervor, indem »® durch —w?® ersetzt wird. 

Sowohl für die quadratische wie für die kubische In- 
terpolation sind — wie bei der linearen Interpolation 
— hier nur zwei Funktionswerte notwendig. Das Her- 
ausziehen des Faktors 1/6 ist wichtig für das Maschinen- 
rechnen zur Einsparung einer Division (falls nicht Ta- 
feln der ‚„‚Everett-Koeffizienten“ zur Hand sind). 

Da manche Funktionen, z. B. 2 

cos, sinx, coshx, sinhx, J,(®), Io(®) 

stets mit ihrer Ableitung zusammen tabelliert sind, er- 
hebt sich die Frage, ob es nicht möglich ist, diesen 
Sachverhalt auszunutzen. Eine solche Formel soll 
nachstehend angegeben werden. 

Dazu gehen wir aus von der Entwicklung 


ul—u 
Ta a 


D) Ö° Yu 
ual—u)(l+u 5. 
TR  ERR; ö° Yıat--- 


1) ö bedeutet die „zentrale Differenz“, also 
Br ( +2) - ( 2). 
(2) = f\ix 5 I 3 


Wir schreiben der Einfachheit halber ö statt ö. 
© 
?®) Hier ist in der üblichen Schreibweise 
h 
= = We (m, + ko). 
») Bsistts=1—u. 


oder, bis auf Glieder vierter und höherer Ordnung 


1— 1 
&+M)=(1—uy+ ung ent 2 vun) 
u(l— u) ee 
Fat 6 Pyın- 


Wegen 


Weeruig. 
ÖYıra = on gryı 
ergibt sich“bei Vernachlässigung der 4. Differenzen 


1 2) 
&n +76 Yan, = 512 + 2) =ÖRyyo- 
Ferner ist s 

r - \ w° rn ZZ 1 
yo -YyltaWy u 1r.-- 


Auch hier können wir den zweiten Summanden als 
klein vernachlässigen und bekommen schließlich 


| sm. 0, 
ot h)=(l—u)y+ u Yıro 9. [y—yı] (1,3). 
Der Fehler beträgt 

1 1 

DTM © -«)  yıai 


er wirdNull für «=, d.h. in diesem Fall haben wir 


kubische Annäherung. Am größten wird der Fehler 
für « = 0,21. Durch Einsetzen finden wir, daß unsere 
Formel bei n-stelliger Genauigkeit brauchbar ist, wenn 
Yyal=l% —3Yyı + 3%— YÜl< 62,510 
Vergleich mit dem bekannten Kriterium für quadra- 
tische Interpolation 
| 6° Yyıa |< 8- 107n 
zeigt den großen Vorteil unserer Formel; man könnte 
sie deshalb eine ‚„semikubische‘‘ nennen. 
Anwendung auf y = sinh & ergibt . 
sinh (2,4%) = (l—u) sinh x, + usinh x, 
w(l— u) 
+0® moon [cosh x, — cosh x]. 


Als Beispiel soll sinh 2,05522 auf 10 Stellen berechnet 
werden. Quadratische Interpolation ist wegen | ö? Yıal 
= 39.107!° nicht mehr möglich. Mit den Tafelwerten 

sinh 2,055 = 3,83937 25185, cosh 2,055 = 3,96746 58541 

sinh 2,056 = 3,84334 19042, cosh 2,056 = 3,97130 67110 
ist dann sinh 2,05522 = 3,84024 54538. 


4 


S a1 1 3 ; 


man beachte ihre Ähnlichkeit mit der linearen Inter-. 
polation! 0 ö 


Als Beispiel soll exp n/5 (r/5 = 0,62831 85307) auf 10 


Stellen berechnet werden. Quadratische Interpolation ist 
wegen | ö° y,,,| = 18.107:° nicht mehr möglich. Mit den 


' Tafelwerten 


exp 0,628 = 1,87385 91108, exp 0,629 = 1,87573 39072 


' erhalten wir exp n/5 = 1,87445 608751). 
? , 


2. Inverse Interpolation. 
Die bekannte Formel fürinverselineare Interpolation 


lautet h 
eu to —! a (21): 


sie setzt voraus, daß man die beiden Funktionswerte yo 
und %, kennt. Dagegen benutzt die von Newton 
angegebene Methode nur einen Funktionswert und die 


‚erste Ableitung: 


ER 2209,02). 


% 


“Diese Gleichung wendet man mit Vorteil auf y = e” 


an. Denn wegen y = y’ erhalten wir aus ihr zur Be- 
rechnung des natürlichen Logarithmus die Formel 


|e=0+ 9.” 1 BD): 


T 


Entsprechend ist für y=e’ 


— = —ye’+1l LA); 
Beim Maschinenrechnen sind diese Formeln nützlich, 
da sie die von Gl. (2,1) geforderte Division einsparen. 
Z.B. bei Benutzung von Hayashis ‚Fünfstelligen 
Tafeln der Kreis- und Hyperbelfunktionen‘‘ kann man 
sie gut gebrauchen; e* und e”” sind dort ja nebenein- 
ander tabelliert. 


Leider findet man in der Literatur nie Formeln für - 


die so oft notwendige inverse quadratische Interpo- 
lation. Es sollen deshalb die Erweiterungen der For- 
meln (2,1) und (2,2) besprochen werden. 

Für die inverse quadratische Interpolation machen 
wir den Ansatz 


2 = % + 6 (Yy— Yo) + Cı (y— Yo) (Yy— Yı)“ 


c, ergibt sich für y = y,. Dagegen bestimmt sich c, mit 
Hilfe eines dritten Funktionswertes; wählen wir für 
diesen Y,, so bekommen wir die Formel 


Il, s MH Yı%ol\ 
BR eg y | Y = 2,5). 


Nimmt man aber y_, als dritten Funktionswert, so er- 
hält unsere Formel für inverse quadratische Interpo- 
lation die Gestalt 


ee 14 YıY Ba -ı 


Yı,% YY-ı 


YıY-ı 


Durch Anwendung beider Formeln bekommt man die 
Fehlerabschätzung; andererseits ergibt ihre Mittelbil- 


*) Man vergleiche diese einfache Rechnung mit der von 
Hayashiim Anhang zu Seinen vielgebrauchten ‚Sieben- 
und mehrstelligen Tafeln‘‘ durchgeführten unter Anwen- 
dung der Newton’'schen Formel für kubische Interpo- 
lation! 


dung einen verbesserten Wert. T ede der beiden Kor 


meln kann auch zur Fehlerberechnung der inversen 
linearen Interpolation benutzt wren. 
- Als a soll die erste Nullstelle der Bessel- 
Funktion J, (x) berechnet werden. Der Tafel entnimmt 
man die Werte ER 
J, (2,403) = 0,0009480931 J, (2,405) =—0,00009 05580 
J, (2,404) =0,00042 86597 J, (2,406) = —0,00060 95598 
Damit ergibt (2,5) bzw. (2,6) 2,40482 55576 5 bzw. 
2,40482 55577 2. Die 10. Stelle ist ungenau, ihr wahrer 
Wert folgt durch Mittelbildung. . 
Nun wollen wir auch die Ableitungen zur inversen 
quadratischen Interpolation heranziehen. Dazu gehen 
wir aus von der Entwicklung 
0 — 


Br? r 
a Kr, 


VI Kerr 


Yı7T7 9 2% %— 


are tert:. 


oder, indem wir die zweiten Ableitungen mit den ersten 
vereinigen, 


en 


a8) 


Setzen wir inden Ableitungen für & den Mittelwert 


1 
FW t ent 5 ein, so folgt mit 


3 
Y%,=(y—y)f (» TE FE: o) ’ 
(2,9) 


? 1 
N, = M-YWf (2 ar zu 0) 
die einfache Formel 


| Yo | 
a Dr 210): 


sie kommt mit nur zwei Funktionswerten aus, hat auch 
die Gestalt der inversen linearen Interpolation, nur 
daß hier die Ableitungen als Korrekturfaktoren auf- 
treten. Die inverse quadratische Interpolation einer 
Funktion ist dadurch auf zwei direkte Interpolationen 
ihrer ersten Ableitung zurückgeführt. 

Die Genauigkeit kann dadurch geprüft werden, daß 
man die Rechnung mit dem interpolierten Wert = & 
noch einmal durchführt. Wir haben also statt (2,9) zu 


nehmen 
$+% 
WW eh 


h,=-m-nf se 


Größere Genauigkeit erreichen wir mit Formel (2,10), 
wenn in den Ableitungen x nicht durch den Mittelwert 


@,11). 


0) 

%+ 3 ersetzt: wird, sondern durch den Wert von x, 
der sich durch inverse lineare Interpolation ergibt; 
dann ist 


Y,= (y—y)f (2 ae alı Sr a) ö | 


0 


Baer (2,12). 
r, = mWf („+3 | 


Bei der Herleitung von (2,10) wurde von dem Glied 
mit der 3. Ableitung noch der halbe Wert erfaßt, so 
daß wir es auch hier wieder mit einer semikubischen 
Interpolation zu tun haben. Besonders für die Um- 
kehrung von Funktionen, die durch Integrale darge- 
stellt werden (wie es oft bei Problemen der Mechanik 
notwendig ist), eignet sich unsere Formel vorzüglich. 


17 und 1,8 hat. 


Auch bei der Auflösung von Gleichungen w 
(2,10) stets weniger Rechenarbeit ha 
5 Beispiele 1 7 NET: per 
I 2 a 
ID HI Wen 


Man findet leicht, daß f (x) = ne Nullstelle zwischen. 
; Ill ergeben (2,9) und aan 


x = 1,763222. Kontrolle mit Ai) ee da 
ein ist. 


- 6. Dezimalstelle um eine Einheit zu 
Als 2. Beispiel sei die seit Fourier oft behandelte 
. Gleichung % 


0 — 22 —5=0 


Wurzel zwischen 2,09 und 2,10 liegt. Damit ergibt 
(2,9) und (2,10) mit einer Division 2,09455 1479; ein 
Wert, dessen Fehler nur 2 Einheiten der 9. Dezimalen 
beträgt. Ye EL ENS 

Aus (2,7) erkennt man, daß für y = x° die inverse 
quadratische Interpolation sofort möglich ist; wir ha- 
ben dann nämlich Se 


m un nn 


‚und damit ist 
er B 
ehem: %(y — 9%) ) (2,13). 


2% (y —Y) + % Yı— 9%) 

Setzen wir in (2,8) die Exponentialfunktion ein, so 
fällt bei den Ableitungen x heraus und es folgt 

Fee ri 

HB. Fr 


w 
rk 
ei. 
also 

=%+ 


ae Er 
Yy—y) + u) en"? | (2,14). 


Damit haben wir eine wunderbar einfache Formel für 
die inverse quadratische Interpolation der e-Funktion 
gefunden; sie unterscheidet sich von der inversen line- 
aren Interpolation ja nur um den für die ganze Tafel 
konstanten Faktor e*”'?. Die Berechnung der natür- 
lichen Logarithmen mit ihrer Hilfe ist weit einfacher, 
als die direkte Interpolation des auf 10 Stellen tabel- 
lierten Logarithmus. 

Als Beispiel soll In 10 auf 10 Stellen genau berechnet 
werden. Einer Tafel der Exponentialfunktion ent- 
nimmt man die Werte 


exp 2,302 = 9,99415 07814 
exp 2,303 = 10,00414 9930 ) 
exp (—0,0005) = 0,99950 01250 


und erhält In 10 = 2,30258 50930. Die Richtigkeit bis 
zur letzten Stelle kann durch neue Rechnung mit dop- 
pelter Intervallbreite, also mit 


exp 2,304 = 10,01415 90846 
exp (—0,001) = 0,99900 04998 


leicht bestätigt werden. 
F.K. Rubbert, 


5) Wegen einer anderen sehr einfachen Methode zur qua- 
dratischen Approximation dritter und beliebiger Wurzeln 
siehe Rubbert, Zur Radizierung mit der Rechen- 
maschine. 


Göttingen. 


Bemerkung zu der Arbeit von L. Cremer [1]: 
Ein neues Verfahren zur Beurteilung der Sta- 
bilität linearer Regelungs-Systeme. 

L. Cremer[l] und A. Leonhard [2] haben 
unabhängig voneinander ein praktisches Verfahren zur 
Prüfung der dynamischen Stabilität eines Systems, 
insbesondere eines Regelkreises angegeben, bei welchem 


.... vektor von A(p) dreht, wenn y von —obis+: 
benutzt. Durch Probieren sei festgestellt, daß eine 7; Kl) pH 


= PW+i 
Oreme: 
nlich eines die 
und Q anzuwendenden en 


Betrachtung des Winkels aus, um den si 
t. a ’ ie 
Veriner hits findet sich beiRouth (l.c.)noch 

ee Stztman 0 
\ p=ıHriy 


worin & und y zwei reelle Veränderliche sind und sieht 
man P und @aie Funktionen won 6 und g ans ma JamIEN NE 
sich in der x,y- zwei Kurven u 

: Pay)=0 


Q(s, = 0 


aufzeichnen; offenbar entspricht jedem Schnittpunkt 


dieser Wurzelkurven [5] eine Wurzel von 


Apd)=Pa,y)+iawy)=0. 


Im Falle der Stabilität Sei die Penn von 
A (p) nurin der tiven bene liegen; die Zweige 
re "Warzelkurvek müssen daher auf der _posi- 


tiven Seite der y-Achse durch die ganze’positive -. 


ebene miteinander abwechseln; auch * Schnitt- 
punkte mit der y-Achse müssen alternieren. Wir er- 
läutern diese Darstellung von Routh durch ein Bei- 
spiel. Es sei 
Ap)=2P®+4m®+3p+1=0. 

Bild1 zeigt den Verlauf der beiden Wurzelkurven 
P(z,y)=22+42 +32 +1—2y(3r+2 = 
Az,y) =y(6r+82+3)—2yP=0 r 


F 


+2 


#1 


OS; 


Bild 1. 
Die Wurzelkurven P(s, W=222+ 42° +320-+1— 2y!3x-- 2)= 0 
aa, Weyer +sc- N —2yV=0 
Ap)=2pP+4p+3p+1=P+iQ=0 
p=s-+iy. 


des Polynoms 


n— 
Die Asymptoten schneiden sich im Punkte — = 


+ 0 —_ 2 
n 3 
und teilen die Ebene in 4n gleiche Winkelräume. 


% 
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Er. 


En; 
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"Man sieht, daß die Schnittpunkte der beiden Wurzel- 


_ kurven in der linken Halbebene liegen und ihre Schnitt- 


punkte mit der y-Achse alterieren. Der durch die zu- 
gehörige Differentialgleichung beschriebene Vorgang 


" ist also stabil. 


Setzt man nun in den Gl. (1) «= 0, so erhält man 
Py...=1—1y=0 a 

am =yß—2m=0.( 9 

. Routh stellt fest, daß diese Gl. (2 Jim Stabilitäts- 

falle reelle Wurzeln haben müssen, die sich gegenseitig 


trennen; d.h. es muß 2 


Y<Y << <H 
sein, wenn %,, %; die beiden Wurzeln der Gl. P (y) = 0 


_ und 0, Ya, Ya die drei Wurzeln der G1.Q (y) = 0 sind. 


Dies ist aber das Kriterium vonCremerundLeon- 
hard. Die Gl. (2) vermitteln die konforme Abbildung 
der y-Achse auf die P, Q-Ebene; in ihr werden die P- 


“ und die Q-Achse von dem Abbild der y-Achse, d.h.. 
der Ortskurve von (A(p)=P+iQ),_, im Sta- 
bilitätsfalle abwechselnd bei den sich trennenden Null- 


stellen von (2) geschnitten; vgl. Bild 2. 


\y: © 


Bild 2. 
Ortskurve A(p)= P(y)tHiQ(y), z=0 
P(W)z-0 = 1—4y' 
aA(Wa=0= v8 —2yY). 


Der Ortsvektor von A (p) dreht sich im positiven 

Sinne um den Winkel 

®=3n, 
wenn y von — » bis + » läuft; an der y (Frequenz)- 
Teilung der Ortskurve liest man die sich trennenden 
„Schnittfrequenzen“ ab. Leonhard, der auf 
Routh nicht Bezug nimmt, weist übrigens über ihn 
hinausgehend darauf hin, daß das auch von ihm mit- 
geteilte Kriterium für die Stabilität auch auf tran- 
szendente charakteristische Gleichungen anwendbar 
bleibt. 

Es mag noch bemerkt werden, daß es für die Be- 
dürfnisse der Regelungstechnik zweckmäßig ist, das 
Problem der Kriterienfüreinen bestimmten 
Lösungstyp der charakteristischen Gleichung von vorn- 
herein nicht auf die Stabilität zu beschränken, sondern 
allgemein nach den Kriterien dafür zu fragen, daß die 
Wurzeln der charakteristischen Gleichung innerhalb 
oder auf dem Rande eines beliebigen, im Endlichen 
linear begrenzten Gebiets der komplexen Ebene liegen, 
also z. B. in der linken Halbebene (Stabilität), in einem 
zur imaginären Achse parallelen Streifen (obere und 
untere Grenze für die Abklingkonstante), in einem 
Streifen parallel zur reellen Achse (Grenzfrequenz), 


 nigenG ra 
gebiet begrenzt, z.B. y=cx für die Abklingdauer, 
für welche das Verhältnis der beiden Wurzelkomp- 


ee 5 


welcher auch anendlieh schmal 


wird (Abklingdauer). Für alle diese Fälle lassen sich 
nach einem einheitlichen Verfahren 


z. B. mit Hilfe des Ca u c h y’schen Integrals[6] für die 
Anzahl der Nullstellen einer analytischen Funktion in . 


einem Gebiet der komplexen Ebene Kriterien angeben, 


deren jedes man geometrisch oder algebraisch form 
lieren kann. Es ist für die Aufstellung dieser Kriterien 


nur erforderlich, in die Polynome (1) an Stelle der im 
Falle des Stabilitätskriteriums zu substituierenden 
‚Gleichung x = 0 für die y-Achse, die Gleichung derje- 

den einzuführen, welche das kritische Wurzel. 


nenten ein Maß ist. h 
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I. Band, 1. Teil: 


Berlin. Hermann Schmidt. 


Anschauliches zur Konvergenz des Iterations- 
verfahrens von Steffensen. 


Hat man die Nullstellen einer Gleichung 


wor) 
zu bestimmen, so konvergiert das übliche Iterations- 
verfahren nicht, wenn in der Umgebung der Nullstelle 
|p’(x)|>1 ist. Für diesen Fallhaben Steffensen!) und 
anscheinend unabhängig vonihm Lemaitre?) ein Ver- 
fahren vorgeschlagen, das darauf hinauskommt, daß 
man die Kurve p(x) durch eine Sekante ersetzt. Man 
iteriert zweimal in der üblichen Weise, bildet also 
2,= pa); = plX,), 
legt durch die beiden Punkte P,[x,, p(x,)] und 
P,[x,, 9(%,)] eine Gerade, die die Gleichung 
% —% 
Y=% Be (e— 2) | 
hat und bringt diese mit der Geraden y=r zum 
Schnitt. Der Schnittpunkt gibt den neuen Näherungs- 
wert 


um, _ An) 
EEE 
Das Verfahren ist anwendbar sowohl wenn |p’(z)|>1 
wie wenn |p’(x)| <1 ist, wenn also auch das gewöhn- 
liche Iterationsverfahren konvergiert. Durch Wieder- 
holung kann man der Nullstelle beliebig nahe kommen. 
Steffensen gibt mittels der allgemeinen Newton- 
schen Interpolationsformel eine Fehlerabschätzung des 
Verfahrens. Deutet man es in der obigen anschaulichen 
Form, so lassen sich aus den Figuren eine Reihe Aus- 
sagen über die Konvergenz machen. Dabei sind 
folgende Fälle zu unterscheiden: FEN 
Fall I: sign (x) = sign 9 («)= +. Ist wie in 
Bild 1a 9’(x)<1, so ist der neue Näherungswert immer 


2) J.F.Steffensen: Skandinavisk Aktuarietidskrift 16 (1938) 


64— 72. i 
2) G.Lemaitre: Akad. Belg. Bull. Cl. Sei. (5) 28 (1942) 


847 —354, 815 —825. 


g dacht werden kann 
(Aperiodizität) oder in einem Winkelraum, der durch 
vom Nullpunkt ausgehende Halbstrahlen begrenzt 


ı > > ’ y; 
R En ten als der 
 vergiert immer, 


man im ersten Fall mit einem Wert 7, <x, und im 
vn‘ zweiten mit einem a >80 Anal der hit 


ai En Das Verfa nre pi: 
wenn man mik einem Wert & 


beginnt. Ist dagegen wie in Bild 1b 0 <p(x) — £ 
muß man mit einem Wert &ı<&o ee 


Bild 1. sign 9’ (x) =signe’(d)= +, 


# 


Ber beiden Geraden soweit herausrücken, daß der neue 
Näherungswert, der auf der anderen Seite der Nullstelle 
liegt, schlechter als der Ausgangswert ist; insbesondere 
kann das leicht der Fall sein, wenn 9 ’(z) sich nur 
wenig von] unterscheidet. Je mehr 9’(x) von 1 ab- 
weicht, um so bessere Konvergenz wird man bekommen. 
Sie wird bei kleinem 9’’(x) besser als bei großem sein. 

Fall II: sign 9 (x) =— sign p’’(x) = +. Im Falle, 
wo 9'(2)>1 ist Bit 2a) ist Fr neue ee 
kleiner als x,; ist dagegen 0 <op’(x) <1,ist er größer 
alsz,.. Die Annäherung erfolgt bei Ww iederholung 
monoton, wenn man im ersten Fall mit 21 <x,, im 
zweiten, wenn man mit x,>x, beginnt. I übrigen 
gilt das zum Fall I Gesagte. 


BETEN ER 2 au nz 2% 
Bild 2. sign p’ (7) = — sign’ (2) = +: 


AR 


Bl nn BER Eu E —nT 
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Bild 3, Bild 4. 
sign 9’ (7) = — sign "(= —. sign 9’ (2) = sign 9’ (= —. 
FallllI: sign p' HR — sign 9" (x) = — (Bild 3). 


Fall IV: sign 9'(x) = sign p’(x) = — (Bild 4). 

In diesen beiden Fällen ist es gleichgültig, ob 
|p’ (#)|<1 oder |p’(x)| > 1 ist und ob man mit einem 
Wert x, <x, oder x,>x, beginnt. Im Falle II ist 


rege 


s der 


TUNGRÖRVERAEEEE =. 


rungs Formänderungen % 
ziert), der hrindende Vol umdehnung beitzt. 
In Anl Be a re Vorgehen bei 2 


matrix Eh wobei sich der Faktor aus einer ku- 
bischen Gleichung als V: 


erzerrungsinvarian x 
Die genannte Konstruktion erfüllt tatsächlich 2 


derung, daß der Deviator bei i 

-änderungen in den gewöhnlichen In Devibiee Ahergk 
Es erscheint mir jedoch wesentlich, darauf hinzuweisen 
daß der Moufangsche Deviator eine andere ebenso 
selbstverständliche us nicht erfüllt, die zu 
stellen ist, wenn man ihn als der Gestalt- 
änderung allein ansehen soll: Wenn zwei ne 
sich nur durch eine Ähnlichkeitsstreckung uniernchölz \ 
den, so müssen die zugehörigen Deviatoren an _ 

Um dies nachzuprüfen, ist es am einfac 


Trefftzschen Verzerrungstensor direkt durch e | 


geometrische Ver matrix darzustellen, wobei wir 
uns der Einfachheit halber, ebenso wie Moufang 
in ihren wesentlichen Ausführungen es tat, auf den Fall 

kartesischer Koordinaten beschränken. 
Die Matrix, die die Deformation der Umgebung eines 
Punktes x — (a1, x?, z°) aufseinen Bildp 2422 

ee ee sei mit D® bezeichnet: 

IK = DI en (1). 


Es ist also 


ds? = dy? = Didy’- Did! = Ar DIDI 


=dy.(BDIar, 


- © wenn der Querstrich über einer Matrix die hierzu an 


der Hauptdiagonale Gespiegelte angibt und Produkte 
in funktioneller Schreibweise von rechts nach links zu 
lesen sind. 


Er Trefftzsche kovariante Verzerrungstensor ist 
also 
1 Er 
S=-ZE-UN).!.... (2). 

Hieraus erhält man sofort die Volumvergrößerung 
dv | —— Den 
ap = |DO= | SD = VIE-2E] ©. 

DerM oufangsche Sue Whatnun as Gestalt 


4=6-,8- Zr.E-I031 ., 


wobei nachMoufangg gemäß (3) aus der Gleich 
|E—2X|= 1, also 3 & “= 


KBSDI+HE|=1..... (5) 


zu bestimmen ist, 


f et Moufang. Z. angew. Math. Mech. 25—27 (1947), 


ARE zum » Moutanaschen Ve AB 
ee - 2 


zus 


RR ng ae 
SS 
2, 


_ Wir gehen nun von 9 zu einer Deformation D, über, 


Ar . die sich von ® durch eine Ähnlichkeitsstreckung A € 
mit A>0 unterscheidet: 27 FR 


2 DI A: D. 


Dann ist 


ee: 
(9,31 = (Od. 


Also wird nach (4), der neue Deviator: 


ri ern LER 
u 
mit einem gemäß: 
IA rAEl- IE RDI+mEI-1M. 


Soll nun A = X, sein, so muß nach (4) und (6) gelten: 
a 
(—pME- 1-2) (DD) 28), 


Bei A=# 1 müßte nach (8) D Dein Vielfaches u? € 
der Einheitsmatrix, d.h. ® selbst eine Ähnlichkeits- 
transformation sein. In diesem Falle wäre nach (4) 
und (5) in der Tat = WU, = 0. 


‚Für eine beliebige Verzerrung A dagegen erfüllt der 
Moufangsche Deviator die. Forderung A= X, 


nicht. Da er demgemäß nicht Repräsentant der er 


staltänderung allein ist, erscheint mir auch seine von 
Moufang vorgeschlagene Verwendung zu einer Ver- 
allgemeinerung des R.Schmi dtschen Verfestigungs- 
gesetzes nicht als gerechtfertigt. 


Das Versagen des Moufangschen Deviators ge- | 


genüber der genannten Inyarianzforderung scheint mir 
noch mehr als seine evtl. Mehrdeutigkeit darauf hinzu- 
weisen, daß er bei der Aufstellung der Verformungs- 
gesetze Schwierigkeiten erzeugen wird. Der haupt- 
sächlichste Grund für seine Ungeeignetheit liegt wohl 
darin, daß man sich bei seiner Aufstell 
das infinitesimale Vorbild gehalten hat, ein Vorwurf, 


der-in gleicher Weise bereits dem Trefftzschen 


Verzerrungstensor gemacht werden kann. 

Es ist jedoch möglich, auch bei beliebigen krumm- 
linigen Koordinaten für den, Verzerrungstensor einen 
Ausdruck anzugeben, bei dem die gewöhnliche Devia- 
torenbildung bereits eine Trennung in Volum- und Ge- 
staltänderung erzeugt, die alle Forderungen erfüllt, und 


bei dessen Benutzung alle Verformungsgesetze tatsäch- 


lich invariant gegenüber beliebigen Koordinatentrans- 
formationen werden. Hierüber soll in einer weiteren 
Arbeit berichtet werden. 
(Die Zuschrift- hat Fräulein Prof. Dr. Moufang 
vorgelegen.) ; 
Haltingen (Kreis Lörrach). H. Richter. 


BUCHBESPRECHUNGEN. 


Dr. F. Reutter, apl. Prof. a. d. Technischen Hoch- 
schule Karlsruhe, Darstellende Geometrie. 
Band I: Orthogonale Zweitafelprojektion. Karlsruhe 
1947. Verlag von G. Braun. 140 S. m. 144 Abb. Preis 
brosch. 7,50 M. : 


Das vorliegende Buch ist nach den Worten des Ver- 
fassers aus Niederschriften von Vorlesungen entstan- 
den, die dieser an der Technischen Hochschule in Karls- 
ruhe gehalten hat. In der Tat wird der Inhalt des Bu- 
ches dem Leser ganz. im Geist einer lebendigen Vor- 
lesung vermittelt, die sich an Hörer der verschiedensten 
Fachrichtungen wendet. Der Verfasser bedient sich 
durchweg in Wort und Bild der anschaulichen Sprache 
des Ingenieurs; langatmige Konstruktionen, Beschrei- 
bungen und Beweisführungen sind ebenso vermieden 
wie Gedankensprünge und logische Lücken. 

Inhaltlich gliedert sich der vorliegende I. Band in 
insgesamt neun Kapitel; die ersten drei Kapitel be- 
schäftigen sich mit den notwendigsten Grundbegriffen 
(Abbildungsmethoden, allgemeine Axonometrie, Affi- 
nität), die übrigen sechs behandeln den Hauptgegen- 
stand des Buches, die senkrechte Zweitafelprojektion 
(Grundgesetze, Vielflache, Zylinder und Kegel, Rota- 
tionsflächen, Flächendurchdringungen, Schrauben- 
linien und -flächen). Durch die Betonung besonderer 


Abschnitte, wie z. B. durch die ausführliche Behand- 


lung der Flächendurchdringungen, ferner auch durch 


die Wahl der Beispiele und Darstellungen selbst bekun- 


det der Verfasser, daß es ihm mit der Anwendbarkeit 


seiner Lehre auch auf praktisch-technische Probleme 
ernst ist. Hierbei bleiben die herangezogenen Objekte 
auch dort, wo aus didaktischen Gründen gestaltlich 


idealisiert oder von technischem Beiwerk abgesehen . 


werden mußte, stets wirklichkeitsnah. Auf einen nur 
äußerlich-technischen Anstrich, wie er vielfach durch 
unmethodisches Anhäufen einseitig technischer Bei- 
spiele erzielt wird, ist bewußt verzichtet. Druck und 
Ausstattung des Buches, das mit zahlreichen und sehr 
gut wiedergegebenen Figuren versehen ist, sind muster- 
gültig. 

Alles in allem ist es dem Verfasser gelungen, in dem 
vorliegenden ersten Band seiner Darst. Geometrie eine 
vorzügliche Einführung in dieses buchmäßig ziemlich 
schwer lehrbare Gebiet zu geben, die sich besonders für 
die Studierenden als Ergänzung der Vorlesungen oder 
zum, Selbststudium eignet, aus der aber auch andere 
Benutzer Gewinn ziehen werden. Möge das Buch, auf 
dessen baldige Fortsetzung wir hoffen, seinen verdien- 
ten Eingang und Anklang finden. 


Darmstadt. H. Graf. 


EINGEGANGENE BÜCHER 


Dr. W. Weizel (o. Prof. an der Uinversität Bonn), 
Einführung indie Physik, Band I: 
Elektrizität und Magnetismus. (2. verb. 
Aufl.). Meyers kleine Handbücher, Bd. 43. 159 S. m. 
155 Abb. Leipzig 1947. Verlag Bibliographisches In- 
stitut. Preis brosch. 2,90 M. 


Dr. Ing. Winfried Oppelt, Grundgesetze der 
Regelung. (Bücher der Technik. Herausgeber: 
Dr. Ing. Alfred Kuhlenkamp). 118 S. m. 32 Bildern 
und 2$Tafeln. Wolfenbüttel-Hannover 1947. Wolfen- 
bütteler Verlagsanstalt G.m.b.H. (Notdruck). Preis 
brosch.-10,— M. 


Dr. W. Lenz (o. ö. Prof. an der Universität Ham- 
burg, Einführungsmathematik für 
Physiker. (Bücher der Mathematik und Natur 
wissenschaften. Herausgegeben von Dr. Henry Poltz.) 


958. m. 45 Bildern. Wolfenbüttel-Hannover 1947.. 
Wolfenbütteler Verlagsanstalt G.m.b.H. (Notdruck). 
Preis brosch. 8,— M. 


Prof. Dr. W. Blaschke, Projektive Geome- 
trie (Bücher der Mathematik und Naturwissenschaf- 
ten. Herausgegeben von Dr. Henry Poltz.) 160 S. m. 
58 Abb. Wolfenbüttel-Hannover 1947. Wolfenbütteler 
Verlagsanstalt G. m. b. H. (Notdruck). Preis brosch. 
13,—M. R 

Prof. Dr. Ing. €. Weber, Festigkeitslehre 
(Bücher der Technik. Herausgegeben von Dr.Ing. Alfred 
Kuhlenkamp). 106 S. m. 204 Abb. Wolfenbüttel-Han- 
nover 1947. Wolfenbütteler Verlagsanstalt G. m. b. H. 
(Notdruck). Preis brosch. 8,80 M. 


(Ausführliche Besprechung eingegangener Bücher 
bleibt vorbehalten.) 


zustarkan - 


Be: 


} 


2 N Ben 
NACHRI 


Friedrich Schilling 80 Jahre. 


Am 9. April 1948 kann der emeritierte o. Professor 
der Mathematik Friedrich Schilling die 
80. Wiederkehr seines Geburtstages begehen, freilich 
fern von Danzig, wo er an der Techn. Hochschule seit 
1904 so segensreich gewirkt hat. Nach furchtbarer 
Flucht Anfang 1945 hat das Ehepaar Schilling schließ- 
lich in Gladbeck/Westfalen bei einer verheirateten 
Tochter ein Unterkommen gefunden, Aber diese 
Schicksalsschläge haben den Jubilarinseiner Schaffens- 
freude nicht gelähmt: demnächst wird wieder eine Mo- 
nographie von ihm erscheinen, reich ausgestattet mit 
schönen, von ihm selbst gezeichneten Figuren. Ihr 
Titel ‚„‚Nichteuklidische Bewegung im hyperbolischen 
Raum“ läßt dieses Buch Schillings als Fortsetzung 
der in der Danziger Zeit entstandenen. Bücher ver- 
wandter Art erkennen, Ihnen allen liegt der gemein- 
same Plan zugrunde, im Rahmen der projektiven Geo- 
metrie die gestaltlichen Verhältnisse der verschiedenen 
Geometrien anschaulich in allen Einzelheiten zu zeigen, 
Diese geometrische Tendenz bekundet der Schüler 
Felix Kleins schon in seiner großen Dissertation 
„Beiträge zur geometrischen Theorie der Schwarz- 
schen s-Funktion‘‘1), auf Grund deren er am 1. März 
1893 in Göttingen die Doktorprüfung bestand. Die von 
Felix Klein so bezeichnete s-Funktion, bekannt 
auch unter dem Namen Schwar zsche Differential- 
invariante, hängt mit der Riemannschen P-Funk- 
tion zusammen; sie bildet die positive Halbebene. auf 
ein von Kreisbogen begrenztes Dreieck ab. Ich kann 


. mich gut erinnern, wie Felix Kleinin seiner faszi- 


nierenden Vorlesung des Wintersemesters 1893/94 über 
die hypergeometrische Funktion ?) wiederholt auf die 
damals im Druck befindliche Dissertation seines frühe- 


ren Assistenten zurückgriff. 
Der Verfasser stand damals im Vorbereitungsdienst 


- für das höhere Schulamt der Provinz Hannover und hat 


sich schnell durch seine glänzende pädagogische, Be- 
geisterung weckende Begabung so ausgezeichnet, daß 
die vorgesetzte Behörde sein Ausscheiden sehr be- 
dıuerte. Schilling war als Assistent für darstellende 
Geometrie an die T.H. Aachen gegangen und hatte 
sieh dort nach Veröffentlichung weiterer an die Disser- 
tation anschließenden Arbeiten habilitiert. Mit seinem 
Ordinarius Friedrich Schur kamer dann nach 
Karlsruhe, wurde dort bald zum a. o. Professor beför- 
dert und zwei Jahre darauf nach Göttingen berufen, 
Bei Eröffnung der T. H. in Danzig 1904 übernahm der 
Jubilar das Ordinariat für Geometrie. In den Jahren 
seiner Göttinger Professur hat dortin besonderer Stärke 
die Neuorganisation des mathematischen Studiums 


‘eingesetzt, wie sie durch die in der Preußischen Prü- 


fungsordnung von 1898 geschaffene Lehrbefähigung 
für angewandte Mathematik gefordert wurde. Hieraus 
sind Schillings Vorträge auf einem Göttinger 
Fortbildungskurs entstanden). Er spricht hier u.a. 
auch von Nomographie, die in jener Zeit in Deutsch- 
land noch recht unbekannt war. Im Anschluß an 
D’Ocagne hatte Schilling in der Göttinger 


ı) Math. Annalen 44 (1894), 8. 161—260. 

®) Im Druck herausgegeben durch O0. Haupt 1938. 
Verlag Julius Springer, Berlin. Vgl. meine ausf. Bespr. 2. 
angew. Math. Mech. 14 (1934) 8. 2541. 

»)InKlein-Riecke über angewandte Mathema- 
tik u. Physik in ihrer Bedeutung für den Unterricht an den 
höh. Schulen. Nebst Erläuterung der diesbezügl. Göttinger 
Universitätseinrichtungen. Leipzig u. Berlin. B. G. Teub- 
ner, 1900. Schillings Gedanken für die Einrichtung von 
Zeichensälen sind heute noch, wo vieles wieder aufgebaut 
werden muß, sehr beachtlich. 
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Er gedenkt dessen dankbar im Vorwort eines seiner b.: 
Bücher. Im Göttinger Mathematischen Studenten- 9 
verein war „derkleineSchilling,derdurchKlein Kr; 
noch groß wird“, sehr beliebt als Mitglied und 


später A.H. 
Welcher Verehrung sich der Jubilar in den en 


Jahren seiner Danziger Tätigkeit erworben hat, = 
die Würdigung erkennen, die sein inzwischen schon Fr 


verstorbener Kollege Sommer ne Jahren dem 
Siebzigjährigen in einer Danziger Zei gewidmet 
hat: ae Sehilling mit a Energie 
geleistet hat für dieAusgestaltung seines Lehrgebietes, 
durch die Schaffung einer rg en Modell- und 
Apparatensammlung für den Unterricht und für 
Forschung, wird noch reiche Früchte tragen und nie 
vergessen werden können“. Aus seiner Danziger Zeit 
sei eine Festrede „Die bildende Kunst und die Geo- 
metrie‘‘ genannt). Von seinen sonstigen vielen Ver- 
öffentlichungen in verschiedenen Zeitschriften, auch 
der Z. angew. Math. Mech., hier zu sprechen, reicht der 
zur Verfügung stehende Raum leider nicht aus, 

„Bei einem Rückblick auf sein Leben darf sich der 
Jubilar wohl sagen, daß es ein Leben war, reich an 
Arbeit, aber auch an Erfolg‘‘, so heißt es in dem Artikel 
Sommers. Wir können das heute erneut aus- 
sprechen und damit die Hoffnung, daßesFriedrich 
Schilling vergönnt ist, noch manche Jahre, ge- 
tragen von der durch Felix Klein 
Tradition, literarisch zu wirken; vielleicht auch schon 
bei der hundertsten Wiederkehr des Gebutrstages von 
Felix Klein am 25. April 1949, 


Frankfurt a.M. W. Lor ey. 
‘)H.A.Schwarz wurde 1892 Nachfolger von 


Weierstraß in Berlin. 
°) Jahresbericht der DMV. 26. 1918, S. 161—176. 


Mit Beginn des Jahres 1948 erscheint unter Leitung 
von L. H. Donnel und St. Timoshenko ein 
neues Referatenblatt für Mechanik „Applied Mechanics 
Reviews“, das sich die gleiche Aufgabe stellt wie das 
Zentralblatt für Mechanik. Es wird herausgegeben 
von „Ihe American Society o£ Mechanical Engineers“ 
mit Unterstüt von „Ihe Office of Naval Research, 
The Engineering Foundation und Illinois Institute of 
Technology‘‘. 


Freiburgi.Br. Dr. H. Bilharz wurde am 
6. 11.1945 zum Dozenten für Mathematik an der 
Universität Freiburg i. Br. ernannt. 


Münsteri.W. Dr.J. Hermes istmit.der Wahr- 
nehmung einer Diätendozentur für Mathematik bei 
auftragt. \ 
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